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ABSTRACT 

We report the discovery of 16 detached M-dwarf eclipsing binaries with J < 16 mag 
and provide a detailed characterisation of three of them, using high-precision infrared 
^| . light curves from the WFCAM Transit Survey (WTS). Such systems provide the most 

' accurate and model-independent method for measuring the fundamental parameters 

£Sl , of these poorly understood yet numerous stars, which currently lack sufficient obser- 

vations to precisely calibrate stellar evolution models. We fully solve for the masses 
and radii of three of the systems, finding orbital periods in the range 1.5 < P < 4.9 
£sq . days, with masses spanning 0.35 — 0.50M Q and radii between 0.38 — O.5OR , with 

uncertainties of ~ 3.5 — 6.4% in mass and ~ 2.7 — 5.5% in radius. Close-companions in 
short-period binaries are expected to be tidally-lockcd into fast rotational velocities, 
resulting in high levels of magnetic activity. This is predicted to inflate their radii by 
K> ' inhibiting convective flow and increasing star spot coverage. The radii of the WTS 

j_j , systems are inflated above model predictions by ~ 3 — 12%, in agreement with the 

Ci ■ observed trend, despite an expected lower systematic contribution from star spots sig- 

nals at infrared wavelengths. We searched for correlation between the orbital period 
and radius inflation by combining our results with all existing M-dwarf radius mea- 
surements of comparable precision, but we found no statistically significant evidence 
for a decrease in radius inflation for longer period, less active systems. Radius inflation 
continues to exists in non-synchronised systems indicating that the problem remains 
even for very low activity M-dwarfs. Resolving this issue is vital not only for under- 
standing the most populous stars in the Universe, but also for characterising their 
planetary companions, which hold the best prospects for finding Earth-like planets in 
the traditional habitable zone. 

Key words: stars: binaries: eclipsing - stars: late-type - stars: fundamental param- 
eters - stars: rotation - stars: magnetic field - surveys. 



1 INTRODUCTION 
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M-dwarfs (M* < 0.6Mq) constitute mor e than seventy per 
cent of the Galactic stellar population (|Henrv et alj|l997h 
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and consequently, they influence a wide-range of astrophys- 
ical phenomena, from the total baryonic content of the 
universe, to the shape of the stellar initial mass function. 
Furthermore, they are fast bec oming a key player in the 



hunt for Earth-like planets (e.g . Nutzman fe Charbonneavj 



l200Sl : lKoppenhoefer et alJl200alLaw et alj|201ll ). The lower 
masses and smaller radii of M-dwarfs mean that an Earth- 
like companion causes a deeper transit and induces a greater 
reflex motion in its host than it would do to a solar analogue, 
making it comparatively easier to detect Earths in the tradi- 
tional habitable zones of cool stars. The inferred properties 
of exoplanet companions, such as their density, atmospheric 
structure and composition, currently depend on a precise 
knowledge of the fundamental properties of the host star, 
such as its mass, radius, luminosity and effective tempera- 
ture at a given age. Yet, to date, no theoretical model of low- 
mass stellar evolution can accur ately reproduce all of the ob- 
served properties of M-dwarf s (|Hillenbrand fe White! 120041 ; 
lLopez-Morales fc Ribasl 120051 ). which leaves their planetary 
companions open to significant mischaracterisation. Indeed, 
the characterisation of the atmosphere of the super-Earth 
around the M-dwarf GJ 1214 seems to depend on the spot 
coverage of the host star l|de Mooii et alj|2012h . 

Detached, double-lined, M-dwarf eclipsing binaries 
(MEBs) provide the most accurate and precise, model- 
independent means of measuring the f undamental proper- 
ties of low-mass stars l|Andersenlll99lf ). and the coevality 
of the component stars, coupled with the assumption that 
they have the same metallicity due to their shared natal 
environment, places stringent observational constraints on 
stellar evolution models. In the best cases, the uncertain- 
ties on th e masses and radii measured using ME Bs can be 
just 0.5% ^Morales et alj|2009l ; iKraus et al.ll201ll ). However, 
since M-dwarfs are intrinsically faint, only a small number 
of MEBs have been characterised so far with suitable ac- 
curacy to calibrate low-mass stellar evolution models, and 
there are even fewer measurements below ~ O.35M0, where 
stellar atmosp heres are thought to trans port energy purely 
by convection dChabrier fc B araffc 1993). 

More worryingly, existing observations show signif- 
icant discrepancies with stellar models. The measured 
radii of M-dwarfs are inflated by 5 — 10% compared 
to model estimates and their eff ective temperatures ap- 
pear too cool by 3 — 5% (see e.g . Lopez-Morales fc Ribasl 

' 20101 ; 



20051; Ribasl 


2004 


Kraus et al.l 


201 lh 



i — 0/0 (see e.g. ILop 
iMorales et "alTl2010l ; 



Torres et al 



This anomaly has been known for 
some time but remains enigmatic. Bizarrely, the two dis- 
crepancies compensate each other in the mass-luminosity 
plane such that current stellar models can accurately re- 
produce the observed mass-luminosity relationship for M- 
dwarfs. Two different physical mechanisms have been sug- 
gested as t he cause of this apparent radius infla tion: i) 



metallicity dBerger et al.ll2006l; lL6pez-Mqralej|2007f) and ii 
magnetic activity jMullan fc MacDonaldlbOOll ; lRibasll200' 



iTorres et al1l200rj ; [chabrier et al.ll2007ri. 

iBerger et al] (|2006l ) and iLopez-Moralesl (|2007f ) used 
interferometrically-measured radii of single, low-mass stars 
to look for correlation between inflation and metallicity. 
Both studies found evidence that inactive, single stars with 
inflated radii corresponded to stars with higher metallic- 
ity, but this did not hold true for active, fast-rotating sin- 
gle stars and further studies could not confirm the result 



l|Demorv et al.l [20091). While metallicity may play a role in 
the scatter of effective temperatures for a given mass (the 
effective temperature depends on the bolometric luminosity 
which is a function of metallicity), it seems unlikely that it 
is the main culprit of radius inflation. 

The magnetic activity hypothesis is steered by the fact 
that the large majority of well-characterised MEBs are in 
short (< 2 day) orbits. Such short period systems found 
in the field (i.e. old systems) are ex pected to b e tidally- 
synchronised with circularised orbits (|Zahnl [l977f h The ef- 
fect of tidal-locking is to increase magnetic activity and is 
a notion that is supported by observations of synchronous, 
rapid rotation rates in MEBs, a majority of circular orbits 
for MEBs, plus X-ray emission and Ha emission from at 
least one of the components. It is hypothesised that in- 
creased magnetic activity affects the radius of the star in 
two ways. Firstly, it can inhibit the convective flow, thus 
the star must inflate and cool to maintain hydrostatic equi- 
librium. IChabrier et al.l (|2007f ) modelled this as a change in 
the convective mixing length, finding that a reduced mix- 
ing length could account for the inflated radii of stars in 
the partially-radiative mass regime, but it had negligible ef- 
fect on the pred icted radii of stars i n the fully-convective 
regime. However, I Jackson et all (|2009T ) showed that the radii 
of young, single, active, fully-convective stars in the open 
cluster NGC 2516 could be inflated by up to 50%, based on 
radii derived using photometrically-measured rotation rates 
and spectroscopically-measured projected rotational veloc- 
ities. This therefore suggests that inhibition of convective 
flow is not the only factor responsible for the radius anomaly. 

The second consequence of increased magnetic activity 
is a higher production of photospheric spots which has a 
two-fold effect: i) a loss of radiative efficiency at the surface, 
causing the star to inflate and ii) a systematic error in light 
curve solutions due to a loss o f circular symmetry c aused by 
a polar distribution of spots. IMorales' et al.1 {201(j ) showed 
that these two effects could account for ~ 3% and — 6% of 
the radius inflation, respectively, with any any remaining ex- 
cess (0 — 4%) produced by inhibition of convective efficiency. 
This however is only under certain generalisations, such as 
a 30% spot coverage fraction and a concentration of the 
spot distribution at the pole. One would perhaps expect the 
systematic error induced by star spots to be wavelength de- 
pendant, such that radius measurements obtained at longer 
wav elength would be c loser to model predictions. 

IKraus et all (|201lh searched for correlation between the 
radius anomaly and the orbital periods of MEBs, to see if the 
data and the models converged at longer periods (~ 3 days) 
where the stellar activity is less aggravated by fast rotation 
speeds. They found tentative evidence to suggest that this 
is the case but it is currently confined to the realm of small 
statistics. Not long after their study, the MEarth project 
uncovered a 41-day, non-synchronised, non-circularised, in- 
active MEB with radius measurements still inflated on av- 
erage by ~ 4%, despite a detailed attempt to account for 
spot-induced systematics l|lrwin et alj|201lh . They suggest 
that either a much larger spot coverage than the 30% they 
assumed is required to explain the inflation, or perhaps that 
the equation of state fo r low-mass stars, desp ite substantial 
progress (see review bv lChabrier et ai]|2005h . is still inade- 
quate. 

Clearly, a large sample of MEBs with a wide-range of 
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orbital periods is key to defining the magnetic activity ef- 
fect and understanding any further underlying physical is- 
sues for modelling the evolution of low-mass single stars. 
This in turn will remove many uncertainties in the prop- 
erties of exoplanets with M-dwarf host stars. With that in 
mind, this paper presents the discovery of many new MEBs 
to emerge from the WFCAM Transit Survey, including a 
full characterisation to reasonable accuracy for three of the 
systems using 4-m class telescopes, despite their relatively 
faint magnitudes (i = 16.7 — 17.6). 

In Section [2] we describe the WFCAM Transit Survey 
(WTS) and its observing strategy, and Section [3] provides 
additional details of the photometric and spectroscopic data 
we used to fully characterise three of the MEBs. In Section[4] 
we outline how we identified the MEBs amongst the large 
catalogue of light curves in the WTS. Sections [5][7] present 
our analysis of all the available follow-up data used to char- 
acterise three of the MEBs including their system effective 
temperatures, metallicities, Ha emission and surface gravi- 
ties, via analysis of low-resolution spectroscopy, their size- 
ratio and orbital elements using multi-colour light curves, 
and their mass ratios using radial velocities obtained with 
intermediate-resolution spectra. These results are combined 
in Section [5] to determine individual masses, radii, effective 
temperatures. We also calculate their space velocities and as- 
sess their membership to the Galactic thick and thin disks. 
Lastly, in Section |5J we discuss our results in the context of 
low-mass stellar evolution m odels and a mass-ra dius-period 
relationship, as suggested bv lKraus et~ai1 (|201ll ). 



2 THE WFCAM TRANSIT SURVEY 

We identified our new MEBs usi ng observations fro m the 
WFCAM Transit Survey (WTS) (|Birkbv et all 1201 lh . The 
WTS in an on-going photometric monitoring campaign that 
operates on the 3.8m United Kingdom Infrared Telescope 
(UKIRT) at Mauna Kea, Hawaii. Its primary and comple- 
mentary science goals are: i) to provide a stringent obser- 
vational constraint on planet formation theories through a 
statistically meaningful measure o f the occurrence ra te of 
hot Jupiters around low-mass stars l|Kovacs et al .1 2012t ) and 
ii) to detect a large sample of eclipsing binaries stars with 
low-mass primaries and characterise them to high enough 
accuracy such that we strongly constrain the stellar evolu- 
tion models describing the planet-hosting M-dwarfs found 
in the survey. The WTS contains ~ 6, 000 early to mid M- 
dwarfs with J ^ 16 mag, covering four regions of the sky 
which span a total of 6 square degrees. 

We combine the large aperture of UKIRT with the 
Wide-Field Camera (WFCAM) infrared imaging array to 
observe in the J-band (1.25/im), near the peak of the spec- 
tral energy distribution (SED) of a cool star. Our observing 
strategy takes advantage of a unique opportunity offered 
by UKIRT, thanks to the highly efficient queue-scheduled 
operational mode of the telescope. Rather than request- 
ing continuous monitoring, we noted there was room for a 
flexible proposal in the queue, which did not require the 
very best observing conditions, unlike most of the on-going 
UKIRT pro grammes that require photometric skies with see- 
ing < 1.3" (|Lawrence et alj|2007h . The WTS is therefore de- 
signed in such a way that there is always at least one target 



field visible and it can observe in mediocre seeing and thin 
cloud cover. We chose four target fields to give us year-round 
visibility, with each field passing within 15 degrees of zenith. 
To select the fields, we co mbined 2MASS photo metry and 
the dust extinction maps of lSchlegel et al.l (1 19981 ') to find re- 
gions of sky that maximised the number of dwarf stars an d 
maximised the ratio of dwarfs to giants l|Cruz et alj|2003h . 
while maintaining E(B — V) < 0.1. We stayed relatively 
close to the galactic plane to increase the number of early 
M-dwarfs, but restricted ourselves to b > 5 degrees to avoid 
the worst effects of overcrowding. 

The survey began on August 05, 2007, and the eclips- 
ing systems presented in this paper are all found in just one 
of the four WTS fields. The field is centred on RA = 19h, 
Dec = +36d, (hereafter, the 19h field), for which the WTS 
has its most extensive coverage, with 1145 data points as of 
June 16, 2011. Note that this field is very close to, but does 
not overlap with, the Kepler field dBatalha et al.|[2006f ). but 
it is promising that recent work showed the giant contam- 
ination in the Kepler field for magnitudes in a comparable 
range to o ur survey was low (7 ± 3% M-giant fraction for 
K P > 14). lMann etaill2012j . 



3 OBSERVATIONS AND DATA REDUCTION 

3.1 UKIRT /WFCAM J-band photometry 

UKIRT and the WFCAM detector provide the survey with 
a large database of infrared light curves in which to search 
for transiting and eclipsing systems. The WFCAM detec- 
tor consists of four 2048 x 2048 18/im pixel HgCdTe Rock- 
well Hawaii-II, non-buttable, infrared arrays that each cover 
13.65' x 13.65' and are separated by 94% of a chip width 
l|Casali et al.ll2007f ). Each WTS field covers 1.5 square de- 
grees of sky, comprising of eight pointings of the WFCAM 
paw print, exposing for a 9-point jitter pattern with 10 sec- 
ond exposures at each position, and tiled to give uniform 
coverage across the field. It takes 15 minutes to observe an 
entire WTS field (9 x 10s x 8+overheads) , resulting in a 
cadence of 4 data points per hour (corresponding to one 
UKIRT Minimum Schedulable Block). Unless there are per- 
sistently bad sky conditions at Mauna Kea, due to our re- 
laxed observing constraints the WTS usually observes only 
at the beginning of the night, just after twilight in > 1" 
seeing when the atmosphere is still cooling and settling. 

The 2-D image processing of the WFCAM observations 
and the generatio n of light curves closely fo llows that of 
llrwin et al.l (|2007h and is explained in detail in lKovacs et al.l 
(| 20121 ). We refer the avid reader to these publications for an 
in-depth discussion of the reduction techniques but briefly 
describe it here. For image processing, we use the automati- 
cally reduced images from the Cambridge Astronomical Sur- 
vey Unit pipeli nqj, which is based on the INT wide-field 
survey pipeline |lrwin fc Lewis! [ioOlh . This provides the 2- 
D instrumental signature removal for infrared arrays includ- 
ing the removal of the dark and reset anomaly, the flat-field 
correction using twilight flats, decurtaining and sky subtrac- 
tion. We then perform astrometric calibration using 2MASS 



1 http:/ /casu. ast.cam.ac.uk/surveys— projects/wfcam/technical 
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stars in the field-of-view, resulting in an astro-metric accu- 
racy of ~ 20 — 50 mas after correcting for field and differ- 
ential distortiorjf]. For photometric calibration, the detector 
magnitude zero-point is derived for each frame using mea- 
surements of stars in the 2MA SS Point Source Cat alogue 
that fall within the same frame i|Hodgkin et alj|2009h . 

In order to generate a master catalogue of source posi- 
tions for each field in the J-band filter, we stack 20 frames 
taken in the best conditions (i.e. seeing, sky brightness and 
transparency) and run our source detection sof tware on the 
stacked image (|lrwinlll985l ; llrwin fc LewislbOOll ). The result- 
ing source positions are used to perform co-located, variable, 
'soft-edged' (i.e. pro-rata flux division for boundary pixels) 
aperture photom etry on all of the time-series images (see 
llrwin et alj|2007h . 

For each of the four WFCAM detector chips, we model 
the flux residuals in each frame as a function of position 
using a 2-D quadratic polynomial, where the residuals are 
measured for each object as the difference between its mag- 
nitude on the frame in question and its median magnitude 
calculated across all frames. By subtracting the model fit, 
this frame-to-frame correction can account for effects such 
as flat-fielding errors, or varying differential atmospheric ex- 
tinction across each fr ame, which can b e significant in wide- 
field imaging (see e.g. llrwin et aLlfeoOTT ). 

Our source detection software flags any objects with 
overlapping isophotes. We used this information in con- 
junction with a morphological image classification flag 
also generated by the pipeline to identify non-stellar or 
blended objects. The plate scale of WFCAM (0.4" /pix) 
is significantly smaller than those of most small aperture, 
ground- based transit survey instrumen ts, such as Super - 
WASP JPollacco et al.ll2006j), H ATNet (|Bakos et all 120041 1 
and TrES ( Dunham et al.ll2004h . and can have the advan- 
tage of reducing the numbers of blended targets, and there- 
fore the numbers of transit mimics, despite observing fainter 
stars. 

The last step in the light curve generation is to per- 
form a correction for residual seeing-correlated effects caused 
by image blending that are not removed by the frame-to- 
frame correction. For each light curve, we model the devi- 
ations from its median flux as a function of the stellar im- 
age FWHM on the corresponding frame, using a quadratic 
polynomial that we then subtract. We note that this method 
addresses the symptoms, but not the cause, of the effects of 
blending. 

Figure [T] shows the per data point photometric preci- 
sion of the final light curves for the stellar sources in the 
I9hr field. The RMS is calculated as a robust estimator us- 
ing as 1.48XMAD, i.e. the equivalent of the Gaussian RMS, 
where the MAD is t he median of the absolute deviations 
jfioaglin et alJll983h . The upturn between J ~ f 2 — f 3 mag 
marks the saturation limit, so for our brightest objects, we 
achieve a per data point precision of ~ 3 — 5 mmag. The blue 
solid line shows the median RMS in bins of 0.2 mag. The 
median RMS at J = f6 mag is ~ 1% (~ 10 mmag), with 
a scatter of ~ 0.8 — f .5%, and only 5% of sources have an 
RMS greater than 15 mmag at this magnitude. Hence, for 
the majority of sources with J ^ 16 mag, the precision is in 
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Figure 1. The RMS scatter per data point of the WTS light 
curves as a function of WFCAM J magnitude, for sources in the 
19hr field with stellar morphological classification. The RMS is a 
robust estimator calculated as 1.48 X the median of the absolute 
deviations. We achieve a per data point photometric precision of 
3 — 5 mmag for the brightest objects, with a median RMS of ~ 1% 
for J = 16 mag. Saturation occurs between ~ 12 — 13 mag as it 
varies across the field and with seeing. The dashed red horizontal 
line at 3 mmag marks the limit of our photometric precision. The 
blue solid curve shows the median RMS in bins of 0.2 mag. The 
red stars show the positions of the 16 WTS 19hr field MEBs. The 
shorter period MEBs sit higher in the plot. RMS values are given 
in Table [ATI 



theory suitable for detecting not only M-dwarf eclipsing bi- 
naries but a lso transits of mid-M stars by planets with radii 
~ lR e (seeH ovacs et al.ll 20121 for the WTS sensitivity to 
Jupiter- and Neptune-sized planets). The 16 new MEBs are 
shown on the plot by the red star symbols. Note that shorter 
period MEBs sit higher on the RMS diagram, but that gen- 
uine longer period MEBs still have RMS values close to the 
median, due to our robust estimator and the long observing 
baseline of the survey. 

For the MEB light curves characterised in this paper, 
we perform an additional processing step, in which we use 
visual examination to clip several clear outlying data points 
at non-consecutive epochs. 

The WTS J-band light curve data for the MEBs re- 
ported in this paper are given in Table[T] We have adopted a 
naming system that uniquely identifies each source handled 
by our data reduction process, and thus we refer to MEBs 
characterised in this paper as: 19b- 2-01387, 19c-3-01405, and 
19e-3-08413. The first number in the naming strategy gives 
the Right Ascension hour the target field. The following let- 
ter corresponds to one of the eight pointings that make up 
the whole WTS target field. The number between the hy- 
phens denotes which of the four WFCAM chips the source 
is detected on and the final 5 digits constitute the source's 
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Name 


HJD 


>Avts 

(mag) 


(mag) 


Amg 
(mag) 


FWHM 6 
(pix) 


x c 
(pix) 


y c 

(pix) 


19b- 2-01387 
19b- 2-01387 


2454317.808241 
2454317.820311 


14.6210 
14.6168 


0.0047 
0.0047 


0.0001 
0.0002 


2.17 
2.37 


321.98 
321.74 


211.07 
210.88 



Table 1. The WTS J-band light curves of 19b-2-01387, 19c-3-01405 and 19e-3-08413. Magnitudes are given in the WFCAM system. 
IHodgkin et al ] j2009h provide conversions for other systems. The errors, rrj, are estimated using a standard noise model, including 
contributions from Poisson noise in the stellar counts, sky noise, readout noise and errors in the sky background estimation. a Correction 
to the fram e magn itude zero point applied in the differential photometry procedure. More negative numbers indicate greater losses. See 
llrwin et al, (2007). 6 Median FWHM of the stellar images on the frame. c x and y pixel coordinates the MEB systems on the image, 
derived using a standard intensity- weighted moments analysis. (This table is published in full in the online journal and is shown partially 
here for guidance regarding its form and content.) 



unique sequence number in our master catalogue of WTS 
sources. 

Some sources in the WTS fields are observed multiple 
times during a full field pointing sequence due to the slight 
overlap in the exposed areas in the tile pattern. 19c-3-01405 
is one such target, receiving two measurements for every full 
field sequence. The median magnitude for 19c-3-0 1405 on 
each chip differs by 32 mmag. IHodgkin et ail (|2009l ) claim a 
photometric calibration error of 1.5% for WFCAM thus the 
median magnitudes have a ~ 2a calibration error. The pho- 
tometric calibration uses 2MASS stars that fall on chip in 
question, so different calibration stars are used for different 
chips and pointing. We combined the light curves from both 
exposures to create a single light curve with 893+898 = 1791 
data points, after first subtracting the median flux from each 
light curve. The combined light curve has the same out-of- 
eclipse RMS, 8 mmag, as the two single light curves. The 
other two MEBs, 19b-2-01387 and 19e-3-08413, have 900 
and 899 data points and an out-of-eclipse RMS of 5 mmag 
and 7 mmag, respectively. 

We also obtained single, deep exposures of each WTS 
field in the WFCAM Z, Y, J, H and K filters (exposure 
times 180,90,90,4 x 90 and 4 x 90 seconds, respectively). 
These are used in conjunction with g, r, i and z photometry 
from SDSS DR7 to create SEDs and derive first estimates 
of the effective temperatures for all sources in the field, as 
described in Section f4. II 



3.2 INT/WFC i-band follow-up photometry 

Photometric follow-up observations to help test and refine 
our light curve models were obtained in the Sloan i-band us- 
ing the Wide Field camera (WFC) on the 2.5m Isaac Newton 
Telescope (INT) at Roque de Los Muchachos, La Palma. 
We opted to use the INT's Sloan i filter rather than the 
RGO I-band filter as i) it has significantly less fringing and, 
ii) unlike the RGO filter, it has a sharp cut-off at ~ 8500 
A and therefore avoids strong, time-variable telluric wa- 
ter vapour absorption lines, w hich could induce systematic s 
in our time-series photometry l|Bailer- Jones fc LammllioOSt ) . 
The observing run, between July 18 - August 01, 2010, was 
part of a wider WTS follow-up campaign to confirm plan- 
etary transit candidates and thus only a few windows were 
available to observe eclipses. Using the WFC in fast mode 
(readout time 28 sec. for lxl binning), we observed a full 
secondary eclipse of 19b-2-01387 and both a full primary 



Name 


HJD 


(mag) 


°" m *INT 

(mag) 


19b- 2-01387 
19b- 2-01387 


2455400.486275 
2455400.487652 


-0.0044 
-0.0049 


-0.0034 
-0.0024 



Table 2. INT i-band follow-up light curves of 19b-2-01387 and 
19c-3-08413. Am, fflT are the differential magnitudes with respect 
to the median of the out-of-eclipse measurements such that the 
out-of-eclipse magnitude is rrti INT = 0. The errors, <7i, are the 
scaled Gaussian equivalents of the median absolute deviation of 
the target from the reference at each epoch i.e. ai ~ 1.48 x MAD. 
(This table is published in full in the online journal and is shown 
partially here for guidance regarding its form and content.) 



and a full secondary eclipse of 19e-3-08413. The observations 
were centred around the expected times of primary and sec- 
ondary eclipse, allowing at least 30 minutes of observation 
either side of ingress and egress to account for any uncer- 
tainty in our predicted eclipse times based on the modelling 
of the WTS light curves. In total, we observed 120 epochs 
for the secondary eclipse of 19b- 2-01387 using 60s exposures, 
and 89 and 82 data points for the primary and secondary 
eclipse of 19e-3-08413, respectively, using 90s exposures. 

We reduced the data using custom built idl routines to 
perform the standard 2-D image processing (i.e. bias sub- 
traction and flat-field division). Low-level fringing was re- 
moved by subtracting a scaled super sky- frame. To create 
the light curves, we performed variable aperture photome- 
try using circular apertures with the idl routine APER. The 
sky background was estimated using a 3cr-clipped median 
on a 30x30 pixel box, rejecting bad pixels. For each MEB, 
we selected sets of 15-20 bright, nearby, non-saturated, non- 
blended reference stars to create a master reference light 
curve. For each reference star, we selected the aperture with 
the smallest out-of-eclipse RMS. We removed the airmass 
dependence by fitting a second order polynomial to the out- 
of-eclipse data. 

The INT i-band light curve data is presented in Table(2] 
The RMS of the out-of-eclipse data for the primary eclipse 
of 19b-2-01387 is 4.4 mmag while the out-of-eclipse RMS 
values for the primary and secondary eclipses of 19e-3-08413 
are 5.7 mmag and 7.1 mmag, respectively. 



© 2012 RAS, MNRAS 000, HWl 



6 J. L. Birkby et al. 



HJD 


^ m ffIAC80 


CTm 9IAC80 




(mag) 


(mag) 


2455052.51020 


-0.0417 


0.0290 


2455052.51113 


-0.0091 


0.0301 



Table 3. IAC80 g-band follow-up light curve of 19e-3-08413. 
Am 9lAC80 are the differential magnitudes with respect to the 
median of the out-of-cclipse measurements such that the out- 
of-cclipse magnitude is mg 1ACS0 = 0. The errors, <j g , are those 
computed by the IRAF.PHOT package. (This table is published in 
full in the online journal and is shown partially here for guidance 
regarding its form and content.) 



3.3 IAC80/CAMELOT 5 -band follow-up 
photometry 

We obtained a single primary eclipse of 19e-3-08413 in the 
Sloan g-band filter using the CAMELOT CCD imager on 
the 80cm IAC80 telescope at the Observatorio del Teide in 
Tenerife. The observations were obtained on the night of 
08 August 2009, during a longer run to primarily follow-up 
WTS planet candidates. Exposure times were 60 seconds 
and were read out with lxl binning of the full detector, 
resulting in a cadence of 71 seconds, making a total of 191 
observations for the night. 

The time-series photometry was ge nerated using the 
VAPHOT package^ (|Deeg fc Dovlejl200ll ). The bias and flat 
field images were processed using standard iraf routines in 
order to calibrate the raw science images. The light curve 
was then generated using VAPHOT, which is a series of 
modified iraf routines that performs aperture photometry; 
these routines find the optimum size aperture that maxi- 
mize the signal-to-noise ratio for each star. The user can 
specify whether to use a variable aperture to account for a 
time-variable point-spread-function (e.g. due to changes in 
the seeing) or to fix it for all images. For this data set, we 
fixed the aperture and used an ensemble of 6 stars with a 
similar magnitude to the target to create a master reference 
light curve. Finally, a second order polynomial was fitted to 
the out-of-eclipse data the target light curve to remove a 
long-term systematic trend. 

The g-band light curve is shown in the bottom left panel 
of Figure [5] and the data are given in Table [3] The out-of- 
eclipse RMS for the target is 26.9 mmag, which is higher 
than the follow-up with the INT, due to the smaller telescope 
diameter. 



3.4 WHT low-resolution spectroscopy 

We carried out low-resolution spectroscopy during a wider 
follow-up campaign of the WTS MEB and planet candidates 
on several nights between July 16 and August 17, 2010, us- 
ing the William Herschel Telescope (WHT) at Roque de Los 
Muchachos, La Palma. These spectra allow the identifica- 
tion of any giant contaminants via gravity sensitive spectral 
features, and provide estimates of the effective system tem- 



3 http:/ /www. iac.es/galcria/hdccg/ 



Name 


Epoch" 


tint 


Instr. 


Grange 


R 


SNR 






to 




(A) 






19b- 2-01387 


394.71 


300 


ISIS 


6000-9200 


1000 


27 


19c-3-01405 


426.53 


900 


ACAM 


3300-9100 


450 


30 


19e-3-08413 


426.54 


900 


ACAM 


3300-9100 


450 


30 



Table 4. Summary of low resolution spectroscopic observations 
at the William Herschel Telescope, La Palma. a JD- 2455000.0. 



peratures, plus approximate metallicities and chromospheric 
activity indicators (see section [5}. 

We used the Intermediate dispersion Spectrograph and 
Imaging System (ISIS) and the Auxiliary-port Camera 
(ACAM) on the WHT to obtain our low-resolution spec- 
tra. In all instances we used a 1.0" slit. We did not use 
the dichroic during the ISIS observations because it can in- 
duce systematics and up to 10% efficiency losses in the red 
arm, which we wanted to avoid given the relative faintness 
of our targets. Wavelength and flux calibrations were per- 
formed using periodic observations of standard lamps and 
spectrophotometric standard stars throughout the nights. 
Table U summarises our low-resolution spectroscopic obser- 
vations. 

The reduction of the low-resolution spectra was per- 
formed with a combination of iraf routines and custom idl 
procedures. In idl, the spectra were trimmed to encompass 
the length of the slit, bias-subtracted and median-filtered 
to remove cosmic rays. The ACAM spectra were also flat- 
fielded. We corrected the flat fields for dispersion effects us- 
ing a pixel-integrated sensitivity function. The iraf.APALL 
routine was used to identify the spectra, subtract the back- 
ground and optimally sum the flux in apertures along the 
trace. For the ISIS spectrum, wavelength and flux calibra- 
tion was performed with the CuNe+CuAr standard lamps 
and INC flux standard SP2032+248. For ACAM, arc frames 
were used to determine the wavelength solution along the slit 
using a fifth order spline function fit with an RMS ~ 0.2A. 
For flux-calibration, we obtained reference spectra of the 
INC flux standard SP2157+261. 

3.5 WHT/ISIS intermediate-resolution 
spectroscopy 

Modelling the individual radial velocities (RVs) of compo- 
nents in a binary system provides their mass ratio and a 
lower limit on their physical separation. Combining this 
information with an inclination angle determined by the 
light curve of an eclipsing system ultimately yields the true 
masses and radii of the stars in the binary. 

We measured the RVs of the components in our MEBs 
using spectra obtained with the intermediate-resolution, 
single-slit spectrograph ISIS mounted on the WHT. We used 
the red arm with the R1200R grating centred on 8500A, giv- 
ing a wavelength coverage of ~ 8100 — 8900A. The slit width 
was chosen to match the approximate seeing at the time of 
observation giving an average spectral resolution R ~ 9300. 

The spectra were processed entirely with iraf, using 
the CCDPROC packages for instrumental signature removal. 
We optimally extracted the spectra for each object on each 
night and performed wavelength and flux calibration using 
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the semi-automatic kpno.doslit package. Wavelength cal- 
ibration was achieved using CuNe arc lamp spectra taken 
after each set of exposures and flux calibration was achieved 
using observations of spectrophotometric standards. 

3.5.1 Radial velocities via cross-correlation 

The region 8700 — 8850A contains a number of relatively 
strong metallic lines present in M-dwarfs and is free of tel- 
luric absorptio n lines making it amenable for M-dwarf RV 
measurements (jlrwin et alj|2009h . We used the iraf imple- 
mentation of the standard 1-D cross-correlation technique, 
FXCOR, to extract the RV measurements for each MEB com- 
ponent using s ynthetic spectra from the MARCS0 spec- 
tral database ( Gustafsson et all [200.8) as templates. The 
templates had plane-parallel model geometry, a tempera- 
ture range from 2800-5500K incremented in 200K steps, so- 
lar metallicity, surface gravity log(gi) = 5.0 and a 2 km/s 
micro-turbulence velocity, which are all consistent with low- 
mass dwarf stars. The best-matching template i.e. the one 
that maximised the cross-correlation strength of the pri- 
mary component for each object, was used to obtain the 
final RVs of the system, although note that the tempera- 
ture of the best-matching cross-correlation template is not 
a reliable estimate of the true effective temperature. The 
saturated near-infrared Ca II triplet lines at 8498, 8542 and 
8662A were masked out during the cross-correlation. A sum- 
mary of our observations and the extracted radial velocities 
are given in table [S] 



4 IDENTIFICATION OF M-DWARF 
ECLIPSING BINARIES 

4.1 The M-dwarf sample 

It is possible to select M-dwarfs in WTS fields using simple 
colour-colour plots such as those shown in Figure [2] which 
were compiled using our deep WFCAM photometry plus 
magnitudes from SPSS DR7, which has a fortuitous over- 
lap with the 19hr field. Ijones et al.l |l994l ) showed that the 
(i — K) colour is a reasonable estimator for the effective tem- 
perature, however the eclipsing nature of the systems we are 
interested in can cause irregularities in the colour indices, 
especially since the WFCAM photometry was taken at dif- 
ferent epochs to each other and the SDSS photometry. For 
example, a system of two equal mass stars in total eclipse 
result is 0.75 mag fainter compared to its out-of-eclipse mag- 
nitude. We made a more robust sample of M-dwarfs by es- 
timating the effective temperature of each source in the 19h 
field via SED fitting of all the available passbands i.e. SDSS 
g, r, i and 2-band plus WFCAM Z, Y, J, H and if-band. By 
rejecting the most outlying magnitudes from the best SED 
fit, one becomes less susceptible to errors from in-eclipse ob- 
servations. Note that the SDSS it-band magnitudes of our 
redder sources are affected by the known red leak in the u 
filter and are hence excluded from the SED fitting process. 

To perform the SED fitting, we first p ut all the ob- 
serve d p hotometry to the V ega system (see Irlewett et al.l 
120061 and lHodgkin et al]|2009l for conversions). Although the 
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Figure 2. Colour-colour plots of the sources in one of the WF- 
CAM pointings for the 19hr field (black +), overlaid with the full 
19hr field sample of detached MEB candidates (blue filled circles 
and red filled squares) . The filled red squares mark the three MEB 
systems characterised in this paper. The orange crosses mark the 
M-dwarf candidate sources in the pointing (see Section 14. ljl . The 
triangles mark the masses for the given colou r index, derived from 
the 1 Gyr solar metallicity isochrone of the Baraffc et al.l lll998T ) 
low-mass stellar evolution models. The arrows mark the maxi- 
mum reddening vector, assuming a distance of 1 kpc. 



WFCAM p hotometry is calibrate d to 1.5 — 2% with respect 
to 2MASS (|Hodakin et al.ll2009l ). the 2MASS photometry 
also carries its own systematic error, so we assumed an extra 
3% systematic error added in quadrature to the photomet- 
ric errors for each source to account for calibration errors 
between different surveys. We used a simple \ fitting rou- 
tine to compare the data to a set of solar metallicity model 
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Name 


HJD 


Slit 


tint 


SNR 


Phase 


RVi 


RV 2 








(n x sec) 






(Kin/s ) 




19b- 2-01387 


2455395.55200 


1.2 


2 x 1200 


22.7 


0.1422 


-143.2 


8.0 


19b- 2-01387 


2455396.46471 


0.7 


3 x 600 


6.22 


0.7513 


23.7 


-158.0 


19b- 2-01387 


2455407.52383 


1.0 


3 x 900 


14.0 


0.1314 


-137.9 


-4.2 


19b- 2-01387 


2455407.62644 


1.0 


3 x 1200 


8.0 


0.1998 


-155.3 


25.1 


19b- 2-01387 


2455408.38324 


1.0 


3 x 900 


9.1 


0.7049 


14.5 


-157.6 


19b- 2-01387 


2455408.51689 


1.0 


3 x 1200 


12.8 


0.7941 


15.1 


-153.7 


1 QK 9 01 
±y u-z-uioo / 




1 O 






O S700 
U.o / UU 


q a 
-y .o 


1 9 


19b- 2-01387 


2455409.38673 


1.0 


3 x 1200 


14.3 


0.3745 


-128.4 


-4.8 


19c-3-01405 


2455407.43073 


1.0 


1200 + 630 


6.4 


0.2244 


-62.5 


57.0 


19c-3-01405 


2455407.47937 


1.0 


3 x 1200 


5.3 


0.2343 


-57.0 


52.7 


19c-3-01405 


2455407.58012 


1.0 


3 x 1200 


5.3 


0.2547 


-63.8 


54.6 


19c-3-01405 


2455408.46929 


1.0 


3 x 1200 


6.0 


0.4347 


-21.7 


22.0 


19c-3-01405 


2455409.56881 


1.0 


3 x 1200 


6.0 


0.6573 


47.3 


-52.6 


19c-3-01405 


2455409.68190 


0.8 


3 x 1200 


5.1 


0.6802 


42.5 


-64.4 


19c-3-01405 


2455409.47707 


0.8 


3 x 1200 


7.5 


0.6387 


46.4 


-43.6 


19e-3-08413 


2455408.42993 


1.0 


3 x 1200 


7.1 


0.6640 


108.0 


-46.5 


19e-3-08413 


2455408.56307 


1.0 


3 x 1200 


8.7 


0.7435 


113.1 


-58.4 


19o-3-08413 


2455409.43629 


1.0 


3 x 1200 


8.9 


0.2654 


-24.8 


140.9 


19o-3-08413 


2455409.52287 


0.8 


3 x 1200 


7.5 


0.3171 


-27.6 


125.6 


19e- 3-08413 


2455409.61343 


0.8 


3 x 1200 


7.5 


0.3712 


-9.4 


109.1 



Table 5. Summary of intermediate-resolution spectroscopic observations. All observations were centred on 8500A. 



magnitud es at an age o f 1 G yr from the stellar evolution 
models of lBaraffe et al.l (1 19981 ). We linearly interpolated the 
model magnitudes onto a regular grid of 5 K intervals from 
1739 — 6554 K, to enable a more precise location of the \ 
minimum. If the worst fitting data point in the best \ 2 fit 
was more than a 5cr outlier, we excluded that data point 
and re-ran the fitting procedure. This makes the process 
more robust to exposures taken in eclipse. The errors on the 
effective temperatures include the formal la statistical er- 
rors from the \ 2 fit (which are likely to be under-estimated) 
plus an assumed ±100 K systematic uncertainty. This error 
also takes into account the k nown missing opacit y issue in 
the optical bandpasses in the lBaraffe et all (|l99&t ) models. 



Our M-dwarf sample is conservative. It contains any 
source with an SED effective temperature ^ 4209 K, mag- 
nitude J ^ 16 mag and a stellar class morphology flag (as 
determined by the data reduction pipeline). The maximum 
effective temperature corresponds to a radius of O.66R0 at 
the typical field star age of 1 Gyr , accor ding to the stellar 
evolution models of lBaraffe et al.l (| 19981 ). We opted to re- 
strict our MEB search to J 16 mag because the prospects 
for ground-based radial v elocity follow-up are bleak beyond 
J = 16 mag (I ~ 18 mag, lAigrain et al.ll2007| ) if we wish to 
achieve accurate masses and radii that provide useful con- 
straints on stellar evolution models. We found a total of 
2, 705 M-dwarf sources in the 19hr field. 



Table [6] gives the single epoch, deep photometry from 
SDSS and WFC AM, plus the prope r moti ons from the SDSS 
DR7 database (|Munn et all 120041 , |200S| ) for the systems 
characterised in this paper. Their SED-derived system ef- 
fective temperatures, T^sed are given in Table [7] 



Parameter 19b-2-01387 19c-3-01405 19e-3-08413 



aj2000 19:34:15.5 19:36:40.7 19:32:43.2 

<5j2ooo 36:28:27.3 36:42:46.0 36:36:53.5 

H a cosS (" /yr) 0.023 ± 0.003 -0.002 ± 0.004 0.008 ± 0.004 

W C'/yr) 0.032 ±0.003 -0.001 ± 0.004 -0.007 ± 0.004 

g 19.088 ±0.010 20.342 ± 0.024 20.198 ±0.020 

r 17.697 ±0.006 18.901 ± 0.012 18.640 ± 0.009 

i 16.651 ± 0.004 17.634 ± 0.008 17.488 ± 0.005 

z 16.026 ±0.007 16.896 ±0.012 16.847 ± 0.010 

Z 15.593 ±0.005 16.589 ± 0.007 16.156 ±0.006 

Y 15.188 ±0.006 16.432 ±0.011 15.832 ± 0.008 

J 14.721 ±0.004 15.706 ±0.006 15.268 ± 0.005 

H 14.086 ±0.003 15.105 ± 0.006 14.697 ± 0.005 

K 14.414 ± 0.006 14.836 ± 0.007 14.452 ± 0.006 



Table 6. A summary of photometric properties for the three 
MEBs, including our photometrically derived effective tempera- 
tures and spectral types. The proper motions fi a cosS and fig are 
taken from the SDSS DR7 database. SDSS magnitudes g,r,i and 
z are in AB magnitudes, while the WFCAM Z, Y, J, H and K 
magnitudes are given in the Vega system. The errors on the pho- 
tometry are the photon-counting errors and do not include the 
extra 3% systematic error we add in quadrature when perform- 
ing the SED-fitting. Conversio ns of the WFCAM ma gnitudes to 
other systems can be found in lHodgkin et al I J2009I) . Note that 
the WFCAM if -band magnitude for 19b-2-01387 was obtained 
during an eclipse event and does not represent the total system 
magnitude. 



4-1.1 Interstellar reddening 

The photometry for the 19hr field is not dereddened before 
performing the SED fitting. The faint magnitudes of our M- 
dwarf sources implies they are at non-negligible distances 
and that extinction along the line-of-sight may be signif- 
icant. This means that our M-dwarf sample may contain 
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hotter sources than we expect. At J ^ 16 mag, assuming 
no reddening, the WTS is distance-limited to ~ 1 kpc for 
the earliest M-dwarfs (Mj = 6 mag at 1 Gyr for M OV, 
A/* = 0.6M Q , using the models of iBaraffe et alJll998h . We 
investigated the reddening effect in the direction of the 19hr 
field using a model fo r interstellar extinction presented by 
iDrimmel etafl (|2003h . In this model, extinction does not 
have a simple linear dependency on distance but is instead 
a three-dimensional description of the Galaxy, consisting of 
a dust disk, spiral arms as mapped by HII regions, plus 
a local Orion-Cygnus arm segment, where dust parame- 
ters are constrained by COBE/DIRBE far infrared obser- 
vations. Using this model, we calculate that Av = 0.319 
mag {E(B — V) = 0.103 mag) at 1 kpc in the direction of 
th e 19hr field. We used the conversion factors in Table 6 
of ISchlegel et all (| 1998T ) to calculate the absorption in the 
UKIRT and SDSS bandpasses, finding A g — 0.370 mag, 
A K = 0.036 mag, E(r - i) = 0.065, E(i - z) = 0.059 and 
E(J - H) = 0.032. The reddening affect along the line-of- 
sight to the field thus appears to be small. We show this 
maximum reddening vector as an arrow in Figure [5] 

For the most interesting targets in the WTS (EBs or 
planet candidates), we obtain low-resolution spectra to fur- 
ther characterise the systems and check their dwarf-like na- 
ture (see Section [5]). Effective temperatures based on spec- 
tral analysis suffer less from the effects of reddening effects 
because the analysis depends not only on the slope of the 
continuum but also the shape of specific molecular features, 
unlike the SED fitting. Therefore, the SED effective temper- 
atures are only a first estimate and we will later adopt values 
derived by fitting model atmospheres to low-resolution spec- 
tra of our MEBs (see Section |5T4"1) . 



4.2 Eclipse detection 

We made the initial detection of our MEBs during an auto- 
mated search for transiting planets in the WTS light curves, 
for which we used the B ox-Least-Squares ( BLS) algorithm, 
OCCFIT, as described in lAigrain fc Irwinl (|2004 ). and em- 
ployed by I Miller et al.l (|2008l ). The box represents a periodic 
decrease in the mean flux of the star over a short time scale 
(an upside-down top hat). The in-occultation data points 
in the light curves fall into a single bin, /, while the out- 
of-occultation data points form the ensemble O. This sin- 
gle bin approach may seem simplistic but in the absence 
of significant intrinsic stellar variability, such as star spot 
modulation, it becomes a valid approximation to an eclipse 
and is sufficient for the purpose of detection. Given the rel- 
atively weak signal induced by star spot activity in the J- 
band, we did not filter the light curves for stellar variabil- 
ity before executing the detection algorithm. We ran OC- 
CFIT on the M-dwarf sample light curves in the 19h field. 
Our data invariably suffer from correlated 'red' noise, thus 
we adjust the OCCFIT detection statistic, S, which assesses 
the signifi cance of ou r dete ctions, with the procedure de- 
scribed bv lPont et al.l (|2006T ) to derive a new statistic, S VB a- 
This process is exp lained in detail for OCCFIT detections in 
iMiller et all |2008l V 



4.3 Candidate selection 

To automatically extract the MEB candidates from results 
of running OCCFIT on the M-dwarf sample light curves, we 
required that S la d ^ 5 and that the detected orbital period 
must not be near the common window-function alias at one 
day i.e. 0.99 > P > 1.005 days. This gave 561 light curves 
to eyeball, during which we removed objects with spurious 
eclipse-like features associated with light curves near the 
saturation limit. 

In total, we found 26 sources showing significant eclipse- 
features in the 19h field, of which 16 appear to be detached 
and have full-phase coverage, with well-sampled primary 
and secondary eclipses. The detached MEB candidates are 
marked on the colour-colour plot in Figure [2] by the blue 
filled circles and red filled squares. The orbital periods of 
the MEBs corresponding to the blue filled circles are given 
in Table [All and their folded light curves are shown in Fig- 
ures [XT] and |X2] The MEBs corresponding to the red filled 
squares are the subjects of the remaining detailed analysis 
in this paper. 



5 LOW-RESOLUTION SPECTROSCOPIC 
ANALYSIS 

Low-resolution spectra of our three characterised MEBs, as 
shown in Figure O permit a further analysis of their com- 
posite system properties and provide consistency checks on 
the main-sequence dwarf nature of the systems. 

5.1 Surface Gravity 

ISlesnick et~aH (|2006f ) and lLodieu et ail (|201lh have shown 
that the depths of alkaline absorption lines between 6300 — 
8825A can highlight low surface gravity features in low- 
mass stars. We used the spectral indices Nagigg and Ti07i4o 
to search for any giant star contaminants in the MEBs and 
found that all three MEBs have indices consistent with dwarf 
star gravity. We note that our low-resolution spectra were 
not corrected for telluric absorption, which is prevalent in 
the Nagigg region, and thus our measured indices may not 
be completely reliable. However a visual inspection of the 
spectra also reveals deep, clear absorption by the Nal dou- 
blet at 8183A, 8195A as highlighted in Figure [3] which is 
not seen in giant stars. For comparison, we also observed an 
M4III giant standard star, [R78b] 115, shown at the top of 
Figure [3] with the same set up on the same night. It lacks 
the deep Na I doublet absorption lines found in dwarfs and 
its measured spectral indices are TiO7i40 = 2.0 ± 0.2 and 
Nagigg = 0.97±0.04, which places it in the lo w-surface grav- 
ity reg ion for M4 spectral types in Figure 1 1 of lSlesnick et al.l 
(2006). The gravity-sensitive spectral index values for our 
MEBs are given in Table [7] 

5.2 Metallicity 

The profusion of broad molecular lines in M-dwarf spectra, 
caused by absorbing compounds such a s Titanium Oxide 
and Vanadium Oxide redwards of 6000A (|Kirkpatrick et al.l 
Il99ll ). make it difficult to accurately define the continuum 
level, which complicates M-dwarf metallicity measurements. 
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Figure 3. Low-resolution spectra of our three new MEBs plus 
a known M-giant star (top spectrum) for comparison. The TiO 
absorption band at 7100A signifies the onset of the M-dwarf spec- 
tral types. The dotted vertical lines, from left the right, mark the 
Na I, Ho and the Na I doublet rest wavelengths in air. The Na 
I doublet is strong in dwarfs while the Calcium infrared triplet 
at 8498, 8542 and 8662A is strongest in giants. The deep features 
at 7594 and 7685A are telluric O2 absorption. Ha emission is 
clearly present in all three MEBs. 



However, recent work shows that the relative strengths of 
metal hydride and metal oxide molecular bands in low- 
resolution optical wavelengths can be used to separate 
metal-po or subdwarf s from solar-metallicity systems. For 
example, IWoolf et all |2009l ) provided a set of equal metallic- 
ity contours in the pla ne of the CaH2±C aH3 and Ti05 spec- 
tral indices defined bv lReid et al.l (1 19951 ). and they mapped 
the m etallicity index CTiO/CaH described by iLepine et al.l 
(|2007l ) onto an absolute metallicity scale, calibrated by 
metallicity measurements from well-defined FGK stars with 
M-dwarf c ompanions, albeit with a significant scatter of 
~ 0.3 dex. iDhital et al.l (|201ll ) have refined the coefficients 
for CTiO/CaH after finding a slight bias for higher metallic- 
ity in early M-dwarfs. We measured the CaH2+CaH3 and 
Ti05 indices in our MEB spectra and compared them with 
these works. We found that all three of our systems are con- 
sistent with solar metallicity. The measured values of the 
metallicity-sensitive indices for our MEBs are given in Ta- 
ble El 

One should note that further progress has been 
made in M-dwarf metallicity measurements by moving 
to the i nfrared and using both low-resolution i f -band 
spectra (|Roias-Avala et all |2010| ; iMuirhead et al 



and high-resolution J-band spectra ( Onehag et al.l 



2011 



2011 



iDel Burgo et al.ll201ll ). These regions are relatively free of 
molecular lines, allowing one to isolate atomic lines (such 
as Na I and Ca I) and thus achieve a precise continuum 
placement. However, in the spectra of M-dwarf short pe- 
riod binary systems, one must be aware that the presence 
of double-lines and rotationally-broadened features further 
increase the uncertainty in their metallicity estimates. 



5.3 Ha Emission 

All three of our MEBs show clear Ha emission in their low- 
resolution spectra, although it is not possible to discern if 
both components are in emission. The equivalent widths of 
these lines, which are a measure of the chromospheric ac- 
tivity, are reported in Table [3 where a negative symbol de- 
notes emission. Ha emission can be a sign of youth, but we 
do not see any accompanying low-surface gravity features. 
The strength of the Ha emission seen in our MEBs is com- 
parab le with other close binary systems (e.g. iKraus et al.l 
l201ll ) and thus is most likely caused by high magnetic ac- 
tivity in the systems. None of the systems have equivalent 
widths < — 8A, which places them in the non-active accre- 
tion region of the empirically d e rived accretion criterion of 
iBarrado v Navascues fc Martini (|2003l ). 



5.4 Spectral type and effective temperature 

Our low-resolution spectra permit an independent estimate 
of the spectral types and effective temperatures of the MEBs 
to compare with the SED fitting values. Initially, we assessed 
the spectral types using the Hammef@ spectral-typing tool, 
which estimates MK spectral types by measurin g a set of 
atomic and molecular features (|Covev et al.ll2007h - One can 
visually inspect the automatic fit by eye and adjust the 
fit interactively. For the latest-type stars (K and M), the 
automated characterisation is expected to have an uncer- 
tainty of ~ 2 subclasses. We found that 19b- 2-01387 has a 
visual best-match with an M2V system, while the other two 
ME Bs were visually closest to M 3V systems. M-dwarf stud- 
ies (|Reid et al.lll995l ; lGizislll997h have found that the TiO 5 
spectral index could also be used to estimate spectral types 
to an accuracy of ±0.5 subclasses for stars in the range K7V- 
M6.5V. The value of this index and the associated spectral 
type (SpT) are given for each of our three MEBs in Table [7] 
We find a reasonable agreement between the spectral index 
results, the visual estimates and the SED derived spectral 
types. 

IWoolf fc Wallersteinl (|2006h derived a relationship be- 
tween the CaH2 index and the effective temperatures of 
M-dwarfs in the range 3500K< T eff < 4000 K. Table [7] 
gives the value of this index and the associated effective 
temperatures, labelled T e g (CaH2), for our three MEBs. 
IWoolf fc Wallersteinl |2006l ) do not quote an uncertainty on 
the relationship, so we assumed errors of ±150 K. Within 
the assumed errors, the effective temperatures derived from 
the spectral indices and the SED fitting agree, but the re- 
lationship between the CaH2 index and T c g has not been 
robustly tested for the CaH2 values we have measured. 



http: / / www.astro.cornell.edu / ~kcovey / thchammer.html 
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Instead, we have determined the system effective tem- 
peratures for our MEBs by directly comparing the ob- 
served spectra to cool star model atmospheres using a \ 2 - 
minimisation algorithm. This incorporated the observational 
errors, which were taken from the error spectrum produced 
during the optimal extraction of t he spectra. We used a grid 
of NextGen atmospheric models (jAllard et al.|[l997l ) inter- 
polated to the same resolution as our low-resolution spectra. 
The models had increments of AT c ff = 100K, solar metal- 
licity and a surface gravity log (3) = 5.0 (a typical value 
for early-type field M-dwarfs), and spanned 5000 — 8500A. 
During the fitting, we masked out the strong telluric O2 fea- 
tures at 7594, 7685 A and the Ha emission line at 6563 A as 
these are not present in the models, although we found that 
their inclusion had a negligible affect on the results. All the 
spectra were normalised to their continuum before fitting. 
We fitted the ^-distribution for each MEB with a six-order 
polynomial to locate its minimum. The corresponding best- 
fitting T e fi (atmos., adopted) is given in Table [7] Assuming 
systematic correlation between adjacent pixels in the ob- 
served spectrum, we multiplied the formal la errors from 
the X 2- fit by v3 to obtain the final errors on the system 
effective temperatures. 

From here on, our analysis is performed with system 
effective temperatures derived from model atmosphere fit- 
ting. Although our different methods agree within their er- 
rors, the model atmosphere fitting is more robust against 
reddening effects, even if this effect is expected to be small, 
as discussed earlier. 



6 LIGHT CURVE ANALYSIS 

Light curves of an eclipsing binary provide a wealth of in- 
formation about the system, including its orbital geometry, 
ephemeris, and the relative size and relative radiative prop- 
erties of the stars. We used the eclip sing binary software, 
jktEBOP0 (|Southworth et alJl2004bM . to model the light 
curves of our MEBs. jktEBOP is a mo dified version of 
EBOP (Eclipsing Binary Orbit Progra m; iNelson fc Davis! 
1 19721 ; IPopper fc Etzell Il98ll ; lEtzell ll98Ch . The algorithm is 
only valid for well-detached eclipsing binaries with small 
tidal di stortions, i.e near-sph erical stars with oblateness 
< 0.04 (|Popper fc Etzell Il98lf ). A first pass fit with jktE- 
BOP showed that this criterion is satisfied by all three of 
our MEBs. 

The light curve model of a detached, circularised eclips- 
ing binary is largely independent of its radial velocity model, 
which allowed us to perform light curve modelling and de- 
rive precise orbital periods on which to base our follow-up 
multi-wavelength photometry and radial velocity measure- 
ments. The RV-dependent part of the light curve model is 
the mass ratio, q, which controls the deformation of the 
stars. In our initial analysis to determine precise orbital pe- 
riods, we assumed circular stars, which is reasonable for de- 
tached systems, but the observed mass ratios (see Section]?]) 
were adopted in the final light curve analysis. 

jktEBOP depends on a number of physical param- 
eters. We allowed the following parameters to vary for 

6 http:/ /www. astro. keele.ac.uk/~jkt/ 
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all three systems during the final fitting process: i) the 
sum of the radii as a fraction of their orbital separation, 
(-Ri + R2)/a, where Rj is the stellar radius and a is the 
semi- major axis, ii) the ratio of the radii, k = R2/R1, iii) 
the orbital inclination, i, iv) the central surface brightness 
ratio, J, which is essentially equal to the ratio of the primary 
and secondary eclipse depths, v) a light curve normalisation 
factor, corresponding to the magnitude at quadrature phase, 
vi) ecoscj, where e is the eccentricity and u) is the longitude 
of periastron, vii) esinuj, viii) the orbital period, P and ix) 
the orbital phase zero-point, To, corresponding to the time 
of mid-primary eclipse. The starting values of P and To are 
taken from the original OCCFIT detection (see Section [4.2[) . 
In the final fit, the observed q value is held fixed. The reflec- 
tion coefficients were not fitted, instead they were calculated 
from the geometry of the system. The small effect of gravity 
darkening was determined by fixing the gravity darkening 
coefficients to suita ble values f or stars with convective en- 
velopes (/3 = 0.32) (|Lucvlll967f ). jktEBOP will allow for a 
source of third light in the model, whether it be from a gen- 
uine bound object or from some foreground or background 
contamination, so we initially allowed the third light param- 
eter to vary but found it to be negligible in all cases and thus 
fixed it to zero in the final analysis. 

Our light curves, like many others, are not of suffi- 
cient quality to fit for limb darkening, so we fixed the 
limb darkening coefficients for each component star, jktld 
is a subroutine of jktEBOP that gives appropriate limb 
darkening law coefficients for a given bandpass based on 
a database of coefficients calculated from available stellar 
model a tmospheres. We u sed the PHOENIX model atmo- 
spheres (|Claretl2000l .l 2004]) and the square-root l i mb da rken- 
ing law in all cases. Studies such as Ivan Hammd l|l993h have 
shown that the square-root law is the most accurate at in- 
frared wavelengths. For each star, we assumed surface grav- 
ities of log(<?) = 5, a solar metallicity and micro-turbulence 
of 2 km/s, and used estimated effective temperatures for the 
component stars: [T efr ,i, T eff>2 ] =[3500K, 3450K] for 19b-2- 
01387, [T eff ,i,Teff, 2 ] =[3300K, 3300K] for 19e-3-08413, and 
[T off , 1, T off , 2 ] =[3525K, 3350K] for 19c-3-01405. 

Note that we did not iterate the limb darkening coef- 
ficients with the final derived values of Ti and T2 (see Sec- 
tion [Sj as they only differed by ~ 30 K (< la) from the 
assumed values. This would be computationally intensive to 
do and would result in a negligible effect on the final result. 

The phase-folded J-band light curves for the MEBs and 
their final model fits are shown in Figures |3] [S] and [S] while 
the model values are given in Table [8] 

6.1 Error analysis 

jktEBO P uses a Levenbe rg-Marquardt minimisation al- 
gorithm (|Press et al.l Il992t ) for least-squares optimisation 
of the model parameters; however, the formal uncertain- 
ties from least-squares solutions are notorious for underesti- 
mating the errors when one or more model parameters are 
held fixed, due to the artificial elimination of correlations 
between parameters. jktEBOP provides a method for as- 
sessing the la uncertainties on the measured light curve pa- 
rameters through Monte Carlo (MC) simulations. In these 
simulations, a synthetic light curve is generated using the 
best-fitting light curve model at the phases of the actual ob- 
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Figure 4. 19b- 2-01387 Left top panel: full phase WFCAM J-band light curve. Left bottom panel: the INT/WFC j-band light curve 
at secondary eclipse. The solid red and purple lines show the best-fit from jktEBOP. The blue data points in the smaller panels show 
the residuals after subtracting the model. Right: Parameter correlations from Monte Carlo simulations and histograms of individual 
parameter distributions. The red dashed vertical lines mark the 68.3% confidence interval. There is a strong correlation between the 
light ratio, the radius ratio, and the inclination (which is skewed), indicating the difficulty in constraining the model even with our high 
quality light curves. 
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Figure 5. 19c-3-01405 Left: WFCAM J-band light curve. Lines and panels as in Figure [4] The magnitude scale is differential as we 
have combined light curves from two different WFCAM chips. Right: Monte Carlo results with lines as in Figure [4] Our inability to 
constrain the model with follow-up data results in strong correlation between the radius ratio and light ratio and parameter distributions 
that are significantly skewed. There are also degeneracies in the inclination which is expected given the near identical eclipse depths. 
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Name 


T c ff 


T e ff 


T c ff 


SpT 


Ti05 


CaH2 


CaH3 


TiO 7 i40 


Na 8 i 89 


EW(Ha) 




(SED) 


(atmos., adopted) 


(CaH2) 


(Ti05) 












(A) 


19b-2-01387 


3494 ± 116 


3590 ± 100 


3586 ± 150 


M2.7±0.5 


0.52 


0.52 


0.73 


1.46 


0.89 


-3.2 


19c-3-01405 


3389 ± 110 


3307 ± 130 


3514 ± 150 


M2.8±0.5 


0.50 


0.48 


0.75 


1.60 


0.87 


-4.3 


19e-3-08413 


3349 ± 111 


3456 ± 140 


3569 ± 150 


M2.3±0.5 


0.54 


0.51 


0.73 


1.46 


0.90 


-4.1 



Table 7. A summary of the spectral indices, derived effective temperatures and spectral types (SpT) for the three characterised MEBs. 
The photometric estimates are labelled with (SED). They have the smallest errors, which include the formal uncertainties plus a 100 
K systematic uncertainty, but they potentially suffer from reddening effects and undcr-estimation of the errors. Our ado pted effective 
temp eratures are marked (atmos., adopted). They are derived from comparison with the NextGen model atmosphere spectra jAllard et alj 
1997 ) and are more robu st against reddening effects . The (Ti05) and (CaH2) labels mark values derived from the spectral index relations 
of lReid et all Jl995h and lWooIf fc Wallerstei nll l2006l'l . respectively. We use T c ff (atmos., adopted) for all subsequent analysis in this paper. 
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Figure 6. 19e-3-08413 Left top panel: full phase WFCAM J-band light curve. Left middle panel: INT/WFC i-band light curves of 
a primary and a secondary eclipse. Left bottom panel: IAC80 g-band light curve of a primary eclipse. The solid red, purple, and cyan 
lines show the best-fit from jktEBOP. Right: Parameter correlations from residual permutations, which gave the larger errors on the 
model parameters than the Monte Carlo simulations, indicating time-correlated systematics. There are strong correlations between the 
light ratio, radius ratio and inclination. 



servations. Random Gaussian noise is added to the model 
light curve which is then fitted in the same way as the data. 
This process is repeated many times and distribution of the 
best fits to the synthetic li ght curves provide the l a uncer- 
tainties on each parameter. ISouthworth et al.l |2005l ) showed 
this technique is robust and gives similar results to Markov 



Chain Monte Carlo techniques used by others, under the 
(reasonable) assumption that the best fit to the observa- 
tions is a good fit. jktEBOP can also perform a residual 
permutation (prayer bead) bead error analysis which is use- 
ful for assigning realistic errors in the presence of correlated 
noise (|Southwortbll200g|) . For each MEB, we have performed 
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Parameter 


19b- 2-01387 


19c-3-01405 


19e-3-08413 


WTS J-band 








P (days) 


1.49851768 ± 0.00000041 


4.9390945 ± 0.0000015 


1.67343720 ± 0.00000048 


T (HJD) 


2454332.889802 ± 0.000077 


2454393.80791 ± 0.00022 


2454374.80821 ± 0.00016 


(Ri +R2)/a 


0.17818 ±0.00040 


0.07023 ± 0.00035 


0.1544 ±0.0016 


k 


0.967 ± 0.044 


0.987 ±0.081 


0.782 ± 0.070 


J 


0.9307 ± 0.0043 


0.993 ±0.013 


0.8162 ±0.0084 




88.761 ±0.051 


89.741 ±0.053 


87.59 ±0.26 


e cos d 


—0.00020 ± 0.00017 


0.000060 ± O.OOOO08 


— 0.00014 ± 0.00017 


e sin u) 


-0.0007 ± 0.0026 


-0.0041 ± 0.0059 


0.0112 ± 0.0062 


Normalisation (mag) 


14.64726 ± 0.00017 


0.00003 ± 0.00020 


15.22776 ± 0.00020 


Ri/a 


0.0906 ± 0.0020 


0.0354 ±0.0014 


0.0867 ± 0.0027 


Ri/a 


0.0875 ± 0.0021 


0.0348 ±0.0015 


0.0676 ± 0.0040 




0.871 ± 0.076 


0.97 ±0.15 


0.503 ± 0.090 


e 


0.0066 ± 0.0026 


0.0058 ± 0.0043 


0.0114 ± 0.0062 


U (°) 


268.0 ± 1.7 


180.5 ± 90.9 


91.1 ± 1.2 


a j (mmag) 


5.2 


8.4 


8.7 


INT i-band 








J 


0.8100 




0.63 


o-i (mmag) 


5.7 




12.1 


IAC80 g-band 








J 






0.6455 


<Ti (mmag) 






29.9 



Table 8. Results from the J and i-band light curve analysis. Only perturbed parameters are listed. The light curve parameter errors 
are the 68.3% confidence intervals while the model values are the means of the 68.3% confidence level boundaries, such that the errors 
are symmetric. To corresponds to the epoch of mid-primary eclipse for the first primary eclipse in the J-band light curve. Errors on 
19c-3-08413 are from residual permutation analysis as they were the largest, indicating time-correlated systematics. crj^ give the RMS of 
the residuals to the final solutions, where all parameters in the fit are fixed to the quoted values and the reflection coefficients calculated 
from the system geometry. 



both MC simulations (using 10, 000 steps) and a prayer bead 
analysis. The reported errors are those from the method that 
gave the largest uncertainties. The correlations between the 
parameter distributions from the MC and prayer bead anal- 
ysis are shown in Figures |4j [5] and [6] along with histograms 
of the distributions of individual parameters. The distribu- 
tions are not perfectly Gaussian and result in asymmetric 
errors for the 68.3% confidence interval about the median. 
As we wish to propagate these errors into the calculation 
of absolute dimension, we have symmetrized the errors by 
adopting the mean of the 68.3% boundaries (the 15.85% and 
84.15% confidence limits) as the parameter value and quot- 
ing the 68.3% confidence interval as the ±lcr errors. These 
errors are given in Table [8] for each MEB. 



Our follow-up g- and i-band light curves (where avail- 
able) were used to check the J-band solution by modelling 
them with the derived J-band parameters, but allowing the 
surface brightness ratio and the light curve normalisation 
factor to vary. The limb darkening coefficients were changed 
to those appropriate for the respective g- and i-band and 
the reflection coefficients were again determined by the sys- 
tem geometry. The RMS values of the these fits are given in 
Table [8] along with the derived g- and j-band surface bright- 
ness ratio for completeness. The g- and i-band phase-folded 
data is shown overlaid with the models in Figure [4] and [6] 
We find that the J-band solutions are in good agreement 
with the g— and i-band data. 



6.2 Light ratios 

All three of our MEBs exhibit near equal-depth eclipses, 
implying that the systems have components with similar 
mass. This is promising because it suggest relatively large 
reflex motions that will appear as well-separated peaks in a 
cross-correlation function from which we derive RVs. How- 
ever, it is well-known for systems with equal size components 
that the ratio of the radii, which depends on the depth of 
the eclipses, i s very poorly determined by the light curve 
(|Popper|[l984l ). even with the high photo metric precision 
and la r ge number epochs in the WTS (see lAndersen et al] 
l|l980t ); ISouthworth et~all (|2007aT l for other excellent exam- 
ples of this phenomenon). Conversely, (Ri ± i?2)/a, is of- 
ten very well-constrained because it depends mainly on the 
duration of the eclipses and the orbital inclination of the 
system. The reason that the ratio of the radii is so poorly 
constrained stems from the fact that quite different values 
of R2/R1 result in very similar eclipse shapes. 

Unfortunately, we found that all three of our MEBs pre- 
sented problems associated with poorly constrained R2 /Ri , 
revealed in the initial modelling as either a large skew in 
the errors on the best-fit parameters or best-fit solutions 
that were physically implausible. For example, for 19b-2- 
01387, the initial best-fit gave L2/L1 > 1 and R2/R1 > 1 
while T2/T1 < 1. We know from our low-resolution spec- 
troscopy that these stars are very likely to be ordinary main- 
sequence M-dwarfs and while their exact radii may be under- 
estimated by models, they generally obey the trend that less 
massive stars are less luminous, smaller and cooler. We note 
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that lStassun et all |2007l ) found a temperature reversal in a 
system of two young brown dwarfs where the less massive 
component was hotter but smaller and fainter than its com- 
panion. In their case the more massive component, although 
cooler, had an RV curve and eclipse depth that were con- 
sistent. In our cases, the most massive component (smallest 
if*) comes towards us (blue-shift) after the deepest (pri- 
mary) eclipse, so it must be the more luminous component. 
The uncertainty in our modelling is most likely to due to 
insufficient coverage of the mid-eclipse points. 

However, we can try to use external data as an addi- 
tional constraint in the fit. Some authors employ a spectro- 
scopically derived light ratio as a n independent constraint 
on k in the light curve modelling (|Southworth et al.l l2004al . 
l2007bl ; iNordstrom fc Johansenll 19941 ). jktEBOP allows the 
user to incorporate an input light ratio in the model and 
propagates the errors in a robust way. The input light 
ratio adds a poin t in th e flux array at a specific phase 
jSouthworth et all [jpOTbl V If this is supplied with a very 
small error, the point is essentially fixed. We have tried sev- 
eral methods to estimate the light ratio for each of our three 
systems, although we stress here that none of the estimates 
should be considered as significant. One requires high reso- 
lution spectra to extract precise light ratios, via the analysis 
of the equivalent width ratios of metallic l ines, which will be 
well- separated if observed at quadrature (|Southworth et all 
120051 ). With a high resolution spectrum, one can disentangle 
the components of the eclipsing binary and perf orm spectral 
index analysis on the separate components (e.g. Ilrwin et al.l 
l2007t) . 

19b-2-01387 is our brightest system and subsequently 
has the highest signal-to-noise in our intermediate-resolution 
spectra. The best spectrum is from the first night of observa- 
tions. For this system, we estimated the light ratio in three 
ways: i) by measuring the ratio of the equivalent widths 
of the lines in the Na II doublet (shown in Figure [7|, ii) 
by using the two-dimensional c ross-correlation algorithm, 
TODCOR (jZucker fc Mazehl Il994h . which weights the best- 
matching templates by the light ratio and, iii) by investigat- 
ing the variation in the goodness-of-fit for a range of input 
light ratios in the model. 

For the first method, the iraf.splot task was used to 
measure the equivalent width of the Na I doublet feature 
with rest wavelength 8183. 27A for each star. Note that this 
assumes the components have the same effective temper- 
ature. The ratio was EW{2)/EW{1) = 0.3582/0.4962 = 
0.7219. In the second method, we found that only the spec- 
trum from the first night contained sufficient SNR to en- 
able TODCOR to correctly identify the primary and sec- 
ondary components. It is known that TODCOR does not 
perform as well for systems with similar spectral features 
(ISouthworth fc Clausen! [20071 ) so we do not use it to derive 
RVs for our nearly equal mass systems. The TODCOR esti- 
mated light ratio was L2/L1 — 0.846. In the final method, 
we iterated jktEBOP across a grid of initial light ratios be- 
tween 0.6-1.1, in steps of 0.01, with very small errors while 
allowing all our usual parameters to vary. The resulting \ 2 ~ 
distribution is not well-behaved. There is a local and global 
minimum at L2/L1 — 0.72 and L2/L1 = 0.97, respectively, 
but the global minimum is bracketed on one side by a signif- 
icant jump to a much larger x 2 suggesting numerical issues. 
We opted to use the light ratio derived with TODCOR as the 
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Figure 7. 19b-2-01387: A high signal-to-noisc intermediate res- 
olution spectrum taken near quadrature phase of 19b- 2-01387 in 
the Na II doublet wavelength region which we used to measure 
equivalents widths to estimate the light ratio. The purple vertical 
lines show the rest frame wavelength of the doublet at A8183.27A, 
A8194.8lA. The red lines mark the doublet for primary object 
and the green lines mark the secondary doublet lines, based on 
the RVs derived in Section [7] 



input to the model. This value lies half-way between the two 
minimums of the \ 2 distribution, so we supplied it with a 
~ 15% error to allow the parameter space to be explored, 
given the uncertainty in our the measurement. Our follow-up 
i-band data of a single secondary eclipse also prefers a light 
ratio less than unity, but the lack of phase coverage does 
not give a well-constrained model. The resulting parameter 
distributions, shown in Figure[4] show strong correlation be- 
tween the light ratio and R2/R1 as expected. The resulting 
la error boundary for the light ratio, which is computed 
from k and J, is in broad agreement with the methods used 
to estimate it. 

For 19e-3-08413, we obtained additional i-band photom- 
etry of a primary and secondary eclipse, plus a further pri- 
mary eclipse in the g-band. Here, we have estimated the 
light ratio by fitting our two datasets in these wavebands 
separately, using appropriate limb darkening coefficients for 
the i- and (/-bands in each case, and allowing all our usual 
parameters to vary. We find best-fit values from the i- and 
g—bands of L2/L1 = 0.29 and L2/L1 = 0.36, respectively. 
This confirms a light ratio less than unity, but as the light 
ratio depends on the surface brightness ratio, which in turn 
is wavelength dependent, we adopted L2/L1 = 0.29 with 
input with an error of ±0.5 in the final fit to the J-band 
data. Note we chose to use the i-band value as it is closer in 
wavelength to the J-band and the light curve was of higher 
quality. 

In the case of 19c-3-01405, we could not derive a light 
ratio from the low SNR spectra, nor do we have follow-up i- 
band photometry (due to time scheduling constraints) . The 
eclipses are virtually identical so we supplied an input light 
ratio of L2/L1 — 1.0 with an error of 50%. Unfortunately, 
the final error distributions for the parameters are still quite 
skewed, as shown in Figure [5] 
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6.3 Star spots 

For 19e-3-08413, we found that the residual permutation 
analysis gave larger errors, indicating time-correlated sys- 
tematics. We have not allowed for spot modulation in our 
light curve model thus the residuals systematics may have a 
stellar origin. As mentioned previously, we expect star spot 
modulation in the J-band to be relatively weak as the SED 
of the spot and the star at these wavelength are more simi- 
lar than at shorter wavelengths. It is difficult to test for the 
presence of spots in the g- and i-band data as we do not have 
suitable coverage out-of-eclipse. We only have full-phase out- 
of-eclipse observations in a single J-bandpass therefore any 
physical spot model will be too degenerate between temper- 
ature and size to be useful. Furthermore, our J-band data 
span nearly four years, yet spot size and location are ex- 
pected to change on much shorter timescales, which leads 
to a change in the amplitude and phase of their sinusoidal 
signatures. Stable star spot signatures over the full duration 
of our observations would be unusual. The WTS observing 
pattern therefore makes it difficult to robustly fit simple si- 
nusoids, as one would need to break the light curve into 
many intervals in order to have time spans where the spots 
did not change significantly (e.g. three week intervals), and 
these would consequently consist of few data points. Never- 
theless, we have attempted to test for spot modulation in a 
very simplistic manner by fitting the residuals of our light 
curve solutions as a function of time (t) with the following 
sinusoid: 



f(t) = do + ai sin(27r(i/a 2 ) + a 3 ), 



(1) 



where the systemic level (ao), amplitude (di), and phase 
(03) were allowed to vary in the search for the best-fit, while 
the period (0,2) was held fixed at the orbital period as we 
expect these systems to be synchronised (see Table [11] for 
the theoretical synchronisation timescales). Once the best-fit 
was found, the values were used as starting parameters for 
the idl routine mpfitfun, to refine the fit and calculate the 
errors on each parameter. Table [5] summarises our findings. 

There is evidence to suggest a low-level synchronous si- 
nusoidal modulation in 19b-2-01387 and 19e-3-08413 with 
amplitude ~ 1.8 — 3.5 mmag, but we do not find significant 
modulation for our longest period MEB (19c-3-01405). The 
modulation represents a source of systematic error that if 
modelled and accounted for, could reduce the errors our ra- 
dius measurements. However, with only one passband con- 
taining out-of-eclipse variation, we cannot provide a use- 
ful non-degenerate model. Good-quality out of eclipse mon- 
itoring is required and given that spot modulation evolves, 
contemporaneous observations are needed, preferably taken 
at multiple wavele ngths to constrain the spot temperatures 
|lrwin et al.ll201ll ). It is surprising that the apparent spot 
modulation in our MEBs persists over the long baseline of 
the WTS observations and perhaps an alternate explana- 
tion lies in residual ellipsoid variations from tidal effects be- 
tween the two stars. We note here that our limiting errors 
in comparing these MEBs to the mass-radius relationship in 
Section f9. II are on the masses, not the radii. 



Parameter 


19b-2-01387 


19c-3-01405 


19e- 3-08413 


Ai (km/s J 


90.7 ± 1.6 


55.2 ± 2.2 


iZ.L ± z.U 


K 2 (km/s) 


94.0 ± 2.3 


60.2 ±1.4 


95.2 ± 3.0 


7 (km/s) 


-70.7 ± 1.3 


-4.8 ±2.0 


43.8 ± 1.8 


RMSi (km/s) 


1.8 


3.7 


2.7 


RMS2 (km/s) 


5.4 


2.5 


5.0 


Q 


0.965 ± 0.029 


0.917 ±0.042 


0.757 ±0.032 


asini (Rq) 


5.472 ± 0.083 


11.27 ±0.25 


5.53 ±0.12 


Mi sin 3 i (M ) 


0.498 ±0.019 


0.410 ±0.021 


0.462 ± 0.025 


M 2 sin 3 i (M Q ) 


0.480 ±0.017 


0.376 ±0.023 


0.350 ±0.018 



Table 10. Results from radial velocity analysis. 



7 RADIAL VELOCITY ANALYSIS 

To extract the semi-amplitudes (Ki,K2) and the centre- 
of-mass (systemic) velocity, 7, of each MEB system, we 
modelled the RV data using the idl routine mpfitfun 
l|Markwardtll2009h . which uses the Levenberg-Marquardt 
technique to solve the least-squares problem. The epochs 
and periods were fixed to the photometric solution values 
as these are extremely well-determined from the light curve. 
Circular orbits were assumed (e = 0) for all three systems as 
the eccentricity was negligible in all light curve solutions. We 
fitted the primary RV data first using the following model: 



RVi = 7 - ifi sin(27r^) 



(2) 



where <j> is the phase, calculated from the light curve 
solution, and K is the semi-amplitude. To obtain K2, we 
then fitted the secondary RV data points using the equation 
above, but this time fixed 7 to the value determined from 
the primary RV data. 



RV2 =7 + K 2 sin(2ir(t>) 



(3) 



The errors on each RV measurement are weighted by 
the RV error given by iraf.fxcor and then scaled until the 
reduced \ 2 01 the model fit is unity. The RMS of the resid- 
uals is quoted alongside the derived parameters in Table [lOl 
and is treated as the typical error on each RV data point. 
The RMS ranges from ~ 2 — 5 km/s between the systems 
and for the given magnitudes of our sy stems is the same as 
the predictions of lAigrain et al.l |2007f) who calculated the 
limiting RV accuracy for ISIS on the WHT, when using 1 
hour exposures and an intermediate resolution grating cen- 
tred on 8500A. 

The RV curves for the primary and secondary compo- 
nents of the three MEBs are shown in Figure [8] along with 
the residuals of each fit. The error bars are the scaled errors 
from IRAF.FXCOR and serve as an indicator of the signal-to- 
noise in the individual spectra and the degree of mismatch 
with the best template. 



8 ABSOLUTE DIMENSIONS AND SPACE 
VELOCITIES 

Combining the results of the light curve and RV curve mod- 
elling allows us to derive the absolute masses and radii 
of our MEB components. Table [TT] gives these dimensions 
along with the separations, individual effective tempera- 
tures, surface gravities, and bolometric luminosities for each 
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Name 


Amplitude 


Phase 


7 


Y 2 , 
^v, before 


Y 2 „ 
A v, after 


RMS be£ore 


RMS a f tcI 




(mmag) 




(mmag) 






(mmag) 


(mmag) 


19b- 2-01387 


1.83 ±0.23 


2.53 ±0.012 


0.19 ±0.15 


l.n 


1.04 


5.2 


4.9 


19c-3-01405 


0.22 ±0.27 


-1.5 ±1.3 


0.23 ± 0.20 


0.87 


0.87 


8,1 


8,1 


19e-3-08413 


3.47 ±0.32 


-0.143 ± 0.050 


0.39 ± 0.22 


1.32 


1.19 


7.8 


7.5 



Table 9. Results of modelling the light curve model residuals with the simple sinusoid defined by Equation [TJ to test for the presence 
of spot modulation. The terms 'before' and 'after' refer to the reduced x 2 and RMS values before subtracting the best-fit sine curve and 
then after the subtraction. Note: mmag = 10 -3 mag. The RMS b cforc value for 19e-3-08413 is different to Table [H] as one data point was 
clipped due to it being a significant outlier. 



Parameter 


19b-2-01387 


19c-3-01405 


19e- 3-08413 


Afl (M s ) 


0.498 ± 0.019 


0.410 ±0.023 


0.463 ±0.025 


M 2 (Mq) 


0.481 ± 0.017 


0.376 ±0.024 


0.351 ±0.019 


Ri (R©) 


0.496 ± 0.013 


0.398 ±0.019 


0.480 ± 0.022 


Ri (R©) 


0.479 ± 0.013 


0.393 ±0.019 


0.375 ±0.020 


a (R Q ) 


5.474 ± 0.083 


11.27 ±0.27 


5.54 ±0.12 


log(9i) 


4.745 ± 0.039 


4.851 ±0.055 


4.742 ±0.053 


logfe) 


4.760 ±0.035 


4.825 ±0.064 


4.834 ±0.051 


T eff ,i (K) 


3498 ± 100 


3309 ± 130 


3506 ± 140 


T cff ,2 (K) 


3436 ± 100 


3305 ± 130 


3338 ± 140 


-^bol,l(LQ) 


0.0332 ± 0.0042 


0.0172 ± 0.0031 


0.0314 ± 0.0058 


^bol,2(L0) 


0.0289 ± 0.0037 


0.0166 ±0.0031 


0.0167 ± 0.0033 


Ml )bo] 


8.45 ±0.14 


9.16 ± 0.20 


8.51 ± 0.19 


Af 2 ,boi 


8.60 ±0.14 


9.20 ±0.20 


9.26 ± 0.23 


Virot.sync (km/s) 


16.73 ±0.45 


4.08 ± 0.19 


14.51 ± 0.55 


V2rot,sync (km/s) 


16.73 ± 0.45 


4.01 ± 0.20 


11.31 ± 0.70 


tsync (Myrs) 


0.05 


6.3 


0.1 


tcirc (Myrs) 


2.6 


1480 


4.0 


^adopted (P c ) 


545 ± 29 


645 ± 53 


610 ± 52 


U (km/s) 


-63.6 ±7.0 


-2.4 ±9.0 


30.9 ±8.6 


V (km/s) 


1.0 ±7.8 


1.3 ± 12.2 


-10.2 ± 11.8 


W (km/s) 


-37 ±6.4 


-4.2 ±8.5 


30.1 ±8.1 



Table 11. Derived properties for the three MEBs. Vrot.s ync are 
the rotational velocities assuming the rotation period is synchro- 
nised with the orbital period. t S ync and t c ; rc are the th eoreti- 
cal tidal syn chronisation and circularisation timescales from lZahnl 
dl975l , ll977l) 



binary system. The masses and radii lie within the ranges 
0.35— 0.50Mq and 0.37— O.5R.0 respectively, and span orbital 
periods from 1 — 5 days. The derived errors on the masses 
and radii are ~ 3.5 — 6.4% and ~ 2.7 — 5.5%, respectively. 

Eclipsing binaries are one of the first rungs on the 
Cosmic Distance Ladder and have provided independent 
distance measurements within the local group e.g. to the 
Large Magellanic Cloud and to the Andromeda G alaxy 
(|Guinan et alj|l998l ; liTibas et al.|[2005l ; lBonanodl2007l ). The 
traditional method for measuring distances to eclipsing bi- 
naries is to compute the bolometric magnitude using the 
luminosity, radius and effective temperature found from the 
light curve and RV curve analysis. This is combined with 
a bolometric correction and the system apparent magnitude 
to compute the distance. While this can yield quite accurate 
results, the definitions for effective temperature and the zero 
points for the absolute bolometric mag nitude and the bolo- 
metric correction must be consistent (jBessell et alj 1 19981 ; 
iGirardi et ai1l2002h . However, we have opted to use a dif- 
ferent method to bypass the uncertainties attached to bolo- 
metric corrections. We used jktabsdim (|Southworth et all 



120051 ). a routine that calculates distances using empirical 
relations between surface brightness and effective temper- 
ature. These relations are robustly tested for dwarfs with 
T c ft > 3600 K and there is evidence that they are valid in 
the infrared to ~ 3000 K jKervella et alj|2004l ). The scat- 
ter around the calibration of the relations in the infrared 
is on the 1% level. The effective temperature scales used 
for the EB analysis and the calibration of its relation with 
surface brightness should be the same to avoid systematic 
errors but this is a more relaxed constraint than require d 
by bolometric correction methods (|Southworth et al]|2005l '). 
The infrared J, H and A'-bands are relatively unaffected by 
interstellar reddening but we have shown in Section 14.1.11 
that we expect a small amount. In the distance determi- 
nation, we have calculated the distances at zero reddening 
and at the maximum reddening (E(B — V) = 0.103 at 1 
kpc for early M-dwarfs with J ^ 16 mag). Our adopted dis- 
tance, dadoptcd, reported in Table fTTI is the mid-point of the 
minimum and maximum distance values at the boundaries 
of their the individual errors, which includes the propaga- 
tion of the effective temperature uncertainties. The MEBs 
lie between ~ 550 — 650 pc. 

With a full arsenal of kinematic information (dis- 
tance, systemic velocities, proper motions and positions) 
we can now derive the true space motions, UVW, for the 
MEBs and determine whether they belong to the Galac- 
tic !_dJsk_oj^lialo_st : el^ We used the method 
of I Johnson fc Soderbloml (|l987h to determine UVW values 
with respect to the Sun (heliocentric) but we adopt a left- 
handed coordinate system to be consistent with the litera- 
ture, that is, U is positive away from the Galactic centre, 
V is positive in the direction of Galactic rotation and W is 
positive in the dire ction of the north Gala c tic po le. We use 
the prescription of I Johnson fc Soderblorr] (|l987f ) to propa- 
gate the errors from the observed quantities and the results 
are summarised in Table fTTI 

Figure [9] shows the MEBs in relation to the error ellipse 
for the Galactic young disk as defined by iLeggettl (|l992t ) 
(-20 < U < 50, -30 < V < 0, -25 < W < 10 w.r.t 
the Sun). 19c-3-01405 is consistent within its error with the 
young disk. 19b- 2-01387 i s an o utlier to the young disk cri- 
terion. Instead, ILeggettl (| 19921 ) define objects around the 
edges of the young disk ellipse as members of the young-old 
disk population, which has a sub-solar metallicity (—0.5 < 
[m/H] < 0.0). 19e-3-08413 exceeds the allowed W range 
for th e young disk, despite overlap in the UV plane. ILeggettl 
(|l992h assign these objects also to the young-old disk group. 
This suggest that two of our MEBs could be metal-poor but 
our spectral index measurements in Section [5] are not accu- 
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Figure 9. The UVW space motions with respect to the Sun 
for our MEBs. The e r rors h ave been propagated according to 
Ijohnson fc Soderbloml |i98t]). The s olid ellipses ar e the error el- 
lipses for the young disk defined bv lLegeetti dl992f) . The dashed 
vertical lines in the lower plot mark the W b oundary within which 
the young-old disk population is contained (Lcggctt 1992). 



Figure 8. Primary and secondary RV curves for the MEBs. Top: 
19b-2-01387; Middle: 19c-3-01405; Bottom: 19e-3-08413. 

The solid black circles are RV measurements for the primary star, 
while open circles denote the secondary star RV measurements. 
The solid red lines are the model fits to the primary RVs and 
the dashed green lines are the fits to the secondary RVs, fixed to 
the systemic velocity of their respective primaries. The horizontal 
dotted lines mark the systemic velocities. The error bars are from 
IRAF.FXCOR but are scaled so that the reduced \ 2 °f the model 
fit is unity. They are merely an indication of the signal-to-noise 
of the individual spectra and the mismatch between the template 
and data. Under each RV plot is a panel showing the residuals of 
the best-fits to the primary and secondary RVs. Note the change 
in scale for the y-axis. The typical RV error for each component 
is given in Table ITOl by the RMS of their respective residuals. 



rate enough to confirm this. We would require, for example, 
higher re solution, J-band spe ctra to assess the metallicities 
in detail (|Onehag et al.ll201lf ) . Comparisons with space mo- 
tions of solar neighbourh ood moving groups do not reveal 
any obvious associations |Soderblom fc Mavorll 19931 ). 



9 DISCUSSION 

9.1 The mass-radius diagram 

Figure [TU] shows the positions of our MEBs in the mass- 
radius plane and compares them to literature mass-radius 
measurements derived from EBs with two M-dwarfs, EBs 
with an M-dwarf secondary but hotter primary, eclipsing 
M-dwarf - white dwarf systems, and inactive single stars 
measured by interferometry. We only show values with re- 
ported mass and radius errors comparable to or better than 
our own errors. The solid line marks th e 5 Gy r, solar metal- 
licity isochrone from the iBaraffe et all (|199Sh models (solid 
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line), with a convective mixing length equal to the scale 
height (L m i x = Hp), while the dash-dot line shows the cor- 
responding 1 Gyr isochrone. 

It is clear that some MEBs, both in the WTS and in 
the literature, have an excess in radius above the model pre- 
dictions, and although there is no evidence to say that all 
MEBs disagree with the models, the scatter in radius at 
a given mass is cl ear, indicating a res idual dependency on 
other parameters. iKnigge et al.l l)201ll ) measured the aver- 
age fractio nal radius excess belo w 0.7 Mp), but based on th e 
findings of lChabrier et alj (120071 ) and lMorales et al.l (|2010l ). 
split the sample at the fully-convective boundary to inves- 
tigate the effect of inhibited convection. The dashed line 
in Figure [TO] marks the average radius inflation they found 
with respect to the 5 Gyr isochrone for the fully-convective 
mass region below 0.35Mq and in the partially-convective 
region above (7.9% for > 0.35M Q , but only by 4.5% for 
> O.35M ). The WTS MEBs sit systematically above the 5 
Gyr isochrone but appear to have good agreement with the 
average radius inflation for their mass range. It is interesting 
to note that we find similar radius excesses to the literature 
despite using infrared light curves. At these wavelengths, we 
crudely expect lower contamination of the light curves by si- 
nusoidal star spots signals and less loss of circular symmetry, 
on account of the smaller difference between the spectral en- 
ergy distributions of the star and the spots in the J-band. 
If one could eliminate the ~ 3% systematic errors i n MEB 
radii caused by polar star spots ^Morales et alJlioiOT ) by us- 
ing infrared data, yet still see similar excess, this would be 
evidence for a larger effect from magnetic fields (or another 
hidden parameter) than currently thought. Unfortunately, 
the errors on our radii do not allow for a robust claim of 
this nature, but it is an interesting avenue for the field. 

The components of our new MEBs do not seem to con- 
verge towards the standard 5 Gyr isochrone as they ap- 
proach the fully-convective region. In fact, our lowest mass 
star, which has a mass error bar that straddles the fully- 
convective boundary, is the most inflated of the six com- 
ponents we have measured. The lower panel of Figure [10] 
illustrates this inflation more clearly by showing the radius 
anomaly Robs/ '-Rmodci as a function of mass, as computed 
with the standard 5 Gyr isochrone. The errors on the ra- 
dius anomaly include the observed error on the radius and 
the observed error on the mass (which propagates into the 
value of 7i mo dei), added in quadrature. The spread in radii 
at a given mass is clearer here, and we discuss why stars 
of the same mass could be inflated by different amounts in 
Section f9. 31 by considering their rotational velocities. 

A comparison of the measured radii of all known MEBs 
to the model isochrones shown in Figure \W\ might lead one 
to invoke young ages for most of the systems, because stars 
with < 0.7Mq are still contracting onto the pre-main 
sequence at an age < 200 Myr and therefore have larger 
radii. While you ng stars exist in t he so lar neighbourhood 
(as shown by e.g. Ijeffries Sz Jewelll |l993l ) who found an up- 
per limit of 10-15 young stars within 25pc), it is highly un- 
likely that all of the known MEBs are young. Indeed, the 
derived surface gravities for our MEBs are consistent with 
older main-sequence stars. We see emission of Ha in all three 
systems, which can be an indicator of youth, but close bi- 
nary systems are known to exhibit significantly more activity 
than wide binaries or single stars of the same spectral type 



(see e.g. IShkolnik et al.ll201ol ). We therefore do not have in- 
dependent evidence to strongly associate the inflated radii 
of our MEBs with young ages. 

9.2 The mass-T c ff diagram 

As discussed in Section [Tl there is some evidence for a radius- 
met a llicity correlation jBerger et af] 120061 ; lLopez-MoralesI 
120071 ) amongst M-dwarfs. Model values for effective temper- 
atures depend on model bolometric luminosities, which are 
a function of metallicity. Metal-poor stars are less opaque so 
model luminosities and effective temperatures increase while 
the m odel radii shrink by a small amount (|Baraffe et al.l 
Il998l ). Figure QT] shows our MEBs in the mass-T c ff plane 
plus the same literature systems from Figure \W\ where ef- 
fective temperatures are availab le. The two li n es sh ow the 
standard 5 Gyr isochrone of the iBarafie et al.l l|l998l ) mod- 
els for solar metallicity stars (solid line) and for metal-poor 
stars (dot-dash line). 

The large errors in the mass-T c fi plane for M-dwarfs 
mean that it is not well-constrained. Section [S] has already 
highlighted some of difficulties in constraining effective tem- 
peratures and metallicities for M-dwarfs, but one should 
also note that effective temperatures reported in the liter- 
ature are determined using a variety of different methods, 
e.g. broad-band colour indices, spectral indices, or model at- 
mosphere fitting using several competing radiative transfer 
codes. It also involv es a number of diff erent spectral type - 
T c fi relations, and as lRevle et al.l (|201ll ) have demonstrated, 
these can differ by up to 500 K for a given M-dwarf subclass. 

While the intrinsic scatter in the effective tempera- 
tures at a given mass may be caused by metallicity ef- 
fects, the overall trend is that models predict temperatures 
that are too hot compared to observed values, especially 
below 0.45Mq. Our new MEBs, which we determined to 
have metallicities consistent with the Sun, also conform to 
this trend. Furthermore, several studies of the inflated CM 
Dra s ystem have found it to be metal-poor (|Viti et al.lll997l . 
I2002D . whereas models would suggest it was metal-rich for 
its mass, based on its cooler temperature and larger radius 
(see Table lA"3l for data). In this case, the very precisely mea- 
sured inflated radius of CM Dra cannot be explained by a 
high metallicity effect. In fact, the tentative association of 
two of our new MEBs with the sligh t ly me tal-poor young-old 
disk population defined bv lLeggettl (1 19921 ). would also make 
it difficult to explain their inflated radii using the metallicity 
argument . 

The scatter in the mass-T e ff plane can also arise from 
spot coverage due to the fact that very spotty stars have 
cooler effective temperatures at a given mass, and conse- 
quently larger radii for a fixed luminosity. Large spot cov- 
erage fractions are associated with high magnetic activity, 
which is induced by fast rotational velocities. Table [TT1 gives 
the synchronous rotational velocities of the stars in our 
MEBs along with their theoretical timescales for tidal cir- 
cularisation and synchronisation. Among our new systems, 
19c-3-01405 contains the slowest rotating stars (~ 4 km/s) 
on account of its longer orbital period, and its components 
have stellar radii that are the most consistent with the stan- 
dard 5 Gyr model. The other faster rotating stars in our 
MEBs have radii that deviate from the model by more than 
la. We discuss this tentative trend between radius inflation 
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Figure 10. The mass-radius diagram for low-mass stars. The filled circles show literature MEB values with reported mass errors < 6% 
and radius errors < 6.5%. Also shown are literature values for i) the low-mass secondaries of eclipsing binaries with primary masses 
> 0.6Mq, ii) M-dwarfs found in M-dwarf - white dwarf eclipsing binaries (MD-WD), and iii) radius measurements of si ngle M-dwarfs 
from i nterferometric data. The red squares mark the new WTS MEBs. The diagonal lines show model isochrones from the lBaraffe et all 
1 199811 models tfm/H] = 0, Y = 0.275 and 

-■'mix — Hp) } while the vertical dotted line marks the onset of fully - convective envelope s 
1 Chabrier &c Bar affc 1997). ' The dashed line shows the 5 Gyr isochrone plus the average radius excess found by iKnigge et alj d201 if) , 
assuming a discontinuity at the fully-convective transition. Above O.35M0, the model is inflated by 7.9%, but below it is only inflated by 
4.5%. The bottom panel shows the radius anomaly, Rabs/Rmadel computed using the 5 Gyr isochrone and again the dashed line shows 
the corresponding average radius excess found byKnill. The literature data used in these plots are given in Table lA3l 



and rotational velocity (i.e. orbital period, assuming the sys- 
tems are tidally-locked) in the next section. 



9.3 A mass-radius-period relationship? 

In a recent paper, iKraus et al.l (|201lf ) presented six new 
MEBs with masses between 0.38 — O.59M0 and short orbital 
periods spanning 0.6 — 1.7 days. Their measurements com- 
bined with existing literature revealed that the mean radii of 
stars in systems with orbital periods less than 1 were differ- 
ent at the 2.6a level to those at longer periods. Those with 
orbital periods < 1 day were systematically larger than the 
predicted radii by 4.8 ± 1%, wherea s for p eriods > 1.5 days 
the deviation from the lBaraffe et all (| 19981 ) models are much 
smaller (1.7±0.7%). The implication is that a very short or- 
bital period, i.e. very high level of magnetic activity, leads 



to greater radius inflation, and one then expects the level 
of radius inflation to decrease at longer periods. Figure [12] 
shows the radius anomaly (Rots / Rmodei) as a function of 
period for our new MEBs plus literature values whose re- 
ported errors are compatible with our own measurements 
(<TM obs < 6% and cr.R oba < 6.5 %). We used the 5 G yr, solar 
metallicity isochrone from the iBaraffe et all (|l99Sf ) models, 
with I/ m ix = Hp, to derive the radius anomalies. The models 
were linearly interpolated onto a finer grid with intervals of 
O.OOOIM© , and the model photospheric radii were calculated 

using Rmodei = \J imodoi/47rcrr o 4 ff lncidel . 

Despite the small sample, we have performed an er- 
ror weighted statistical analysis of the period distribution, 
including our new measurem ents, to compare t o the un- 
weighted analysis presented in IKraus et all (|201ll ). Table [T2l 
reports the weighted mean (p.) and weighted sample stan- 
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Figure 11. The mass-T e ff diagram for low mass stars. Two different metallicity isochrones from the lBaraffe et aD dl998f) 1 Gyr models are 
over- plotted to show the effect of decreasing metallicity. The vertical dotted line marks the fully-convective boundary jchabrier fc Baraffe! 
Il997tl . The data used in this plot are given in Table IA3I 



dard deviation (it) of the radius anomaly for three different 
period ranges: i) all periods, ii) periods ^ 1 day and, iii) pe- 
riods > 1 day. The boundary between the 'short' and 'long' 
pe riod samples was c hosen initially to match the analysis 
by iKraus et alj |201 l[ ). A T-test using the weighted mean 
and variances of the short and long period samples shows 
that their mean radii a re distinct po pulations at a 4. Oct sig- 
nificance, in support of lKraus et all 's findings. However, the 
significance level is strongly dependent on the chosen period 
boundary, and is skewed by the cluster of very precisely mea- 
sured values near 1.5 days. For example, a peak significance 
of 4.8(7 is found when dividing the sample at 1.5 days, but 
sharply drops to ~ la for periods of 1.7 days or longer. At 
short periods, it rises gradually towards the peak from la 
at 0.3 days. 

Instead, we have attempted to find a very basic math- 
ematical description for any correlation between radius in- 
flation and orbital period, but we appreciate our efforts are 
hampered by small number statistics. We fitted the distri- 
bution of the radius anomaly as a function of period, using 
first a linear model and then as an exponentially decaying 
function. We used the idl routine mpfitfun to determine an 



Period 


P 




(7 


All 


103.6% 


0.5% 


3.2% 


P sC 1.0 


106.1% 


0.9% 


3.5% 


P > 1.0 


102.6% 


0.4% 


2.4% 



Table 12. A statistical analysis of the mean radius inflation for 
different period ranges, <r is the weighted sample standard devia- 
tion. 



error weighted best-fit and the la errors of the model param- 
eters. The results are reported in Table [131 and the best-fit 
models are over-plotted in Figure 1121 but neither model is 
a good fit (although the exponential fairs moderately bet- 
ter). While there is marginal evidence for greater inflation 
in the shortest period systems, we find that the expected 
convergence towards theoretical radius values for longer pe- 
riod, less active systems is not significantly supported by the 
available observation data. 

There are two pertinent observations worth ad- 
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Figure 12. The radius anomaly as a function of orbital period using the 5 Gyr solar-metallicity isochrone from the lBaraffe et al . (1998) 
models. Our new MEBs are shown by the red open squares. Literature radius anomalies with radius errors < 6.5% are also plotted. The 
errors are a quadrature sum of the measured radius error plus a propagated error from the observed mass which determines the model 
radius. The dashed and dotted lines show the best-fit from a straight-line and exponentially decaying model to the data, respectively. 
The coefficients and goodness of fit for these fits are given in Table [l3l The data used in this plot are given in Table lA3l 



Model a oi a 2 X 2 DOF \l 



^obs/Rmod = a +aiP 1.0401 ±0.0017 -0.000386 ± 0.000086 490.5 46 10.7 

RobsARmod = a + ai e a 2 p 1.0221 ±0.0027 0.089 ± 0.015 -1.57 ±0.33 401.6 45 8.9 



Table 13. Results from an error weighted modelling of the radius anomaly as a function of period, are the coefficients of the models 
and P is the orbital period in days. Neither of these simple models provide a statistically good fit, indicating a more complex relationship 
between the radius anomaly and orbital period. 



dressing, namely the low- mass eclipsing binaries LSPM 
J111 2 +7626 and Kepler- 16 (|lrwin et alj|201lt IPovle et all 
20 111; iBender et al. I l2012h , which were announced after the 



Kraus et alJ l|201ll l study. These systems significantly 



tended the observed orbital period range, with almost iden- 
tical 41-day orbital periods, and both c ontaining one fully- 
conve ctive component (M* ~ 0.35Mn. IChabrier fc Baraffd 
1 19971 ) and one partially convective component (see Ta- 
ble IA3[) . The radius inflation differs significantly between 
these two systems, as shown on the right-hand side of Fig- 
ure [T2j While the more massive, partially-convective com- 



ponent of Kepler-16 is well-de scribed by the 1 Gyr model 
isochrone iBaraffe et al.l (| 19981 ) (see Figure I10[) . the other 
three stars suffer significant radius inflation, with no obvious 
correlation between the amount of inflation and the masses, 
even though one of them is a partially-convective star. This 
residual inflation, particularly for the fully-convective stars 
at long periods, may pose a challenge to the magnetic ac- 
tivity hypothesis as the sole reason for discrepancies be- 
tween models and observations, especially given the ex- 
tremely high-quality measurements of Kepler-16. However, 
one should note that other studies have suggested that the 
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presence of a strong magnetic field can alter the interior 
structure of a low-mass star, such that is pushes the fully- 
convective mass limit for very active stars to low er values 
jMullan fc MacDonaldll200ll ; IChabrier et~al]|200^ , so these 
stars may still suffer from a significant inhibition of convec- 
tive flow. 

The radius anomaly raises concern over the usefulness 
of the known MEBs in calibrating models for the evolu- 
tion of singular M-dwarf stars that are the favoured targets 
of planet-hunting surveys searching for habitable worlds. 
iKraus et alj (|201 lh argue that the high-activity levels in very 
close MEBs make them poor representatives of typical sin- 
gle low-mass stars and that the observed radius discrepan- 
cies should not be taken as an indictment of stellar evolution 
models. However, we have seen that radius inflation remains 
in MEBs systems with low magnetic activity and further- 
more, the inflated components of LSPM Jl 1 12+7626 do not 
exhibit Ha emission that is typically associate d with the 
high a ctivity levels in MEBs with inflated radii. IWest et al.l 
|201ll ) used Ha emission as an activity indicator to deter- 
mine that the fraction of single, active, early M-dwarfs is 
small (< 5%), but increases to 40 - 80% for M4-M9 dwarfs. 
Yet, it may be that the amount of activity needed to in- 
flate radii to the measured values in MEBs is small and 
therefore below the level where observable signatures appear 
in Ha emission. This would then question the reliability of 
Ha emission as an activity indicator, meaning the fraction 
of 'active', single M-dwarfs may be even higher than the 
IWest et all |201ll ) study. Given that these very small stars 
are a ripe hunting-ground for Earth-size planets, we must 
be able to constrain stellar evolution models in the pres- 
ence of magnetic activity if we are to correctly characterise 
planetary companions. We note that even the very precisely- 
calib r ated higher-mass stellar evolution models |Andersenl 
Il99ll ; I Torres et al.l l20ld ) do not reproduce the radii of ac- 
tive stars accurately (see iMorales et al.l (|2009h who found 
4 — 8% inflation in a G7+K7 binary with a 1.3 day orbit). 

In order to establish a stringent constraint on the rela- 
tionship between mass, period and radius, we need further 
measurements of systems that i) include 'active' and 'non- 
active' stars that span the fully-convective and partially- 
convective mass regimes, and ii) a better sampled range of 
orbital periods beyond 5 days to explore systems that are 
not synchronised. We may ultimately find that activity does 
not account for the full extent of t he radius anomaly, and 
as suggested by I Irwin et al.l ( |201lf ) . perhaps the equation 
of state for low-mass stars can still be improved. On the 
other hand, perhaps the importance of tidal effects between 
M-dwarfs in binaries with wider separations has been under- 
estimated, as it has been shown that the orbital evolution 
of M-d warf binary system s is not well-described by current 
models (|Nefs et al.ll2012T ). 



10 CONCLUSIONS 

In this paper, we have presented a catalogue of 16 new low- 
mass, detached eclipsing binaries that were discovered in the 
WFCAM Transit Survey. This is the first time dynamical 
measurements of M-dwarf EBs have been detected and mea- 
sured primarily with infrared data. The survey light curves 
are of high quality, with a per epoch photometric precision 



of 3 — 5 mmag for the brightest targets (J ~ 13 mag), and a 
median RMS of < 1% for J < 16 mag. We have reported the 
characterisation of three of these new systems using follow- 
up spectroscopy from ground-based 2 — 4 m class telescopes. 
The three systems (i — 16.7— 17.6 mag) have orbital periods 
in the range 1.5 — 4.9 days, and span masses 0.35 — O.5OM0 
and radii 0.38 — O.5OR.0, with uncertainties of ~ 3.5 — 6.4% 
in mass and ~ 2.7 — 5.5% in radius. Two of the systems may 
be associated with the young-old disk population as defined 
by iLeggettl (| 1992h but our metallicity estimates from low- 
resolution spectra do not confirm a non-solar metallicity. 

The radii of some of the stars in these new systems are 
significantly inflated above model predictions (~ 3 — 12%). 
We analysed their radius anomalies along with literature 
data as a function of the orbital period (a proxy for activity). 
Our error-weighted statistical analysis revealed marginal ev- 
idence for greater radius inflation in very short orbital pe- 
riods < 1 day, but neither a linear nor exponentially decay 
model produced a significant fit to the data. As a result, we 
found no statistically significant evidence for a correlation 
between the radius anomaly and orbital period, but we are 
limited by the small sample of precise mass and radius mea- 
surements for low-mass stars. However, it is clear that radius 
inflation exists even at longer orbital periods in systems with 
low (or undetectable) levels of magnetic activity. A robust 
calibration of the effect of magnetic fields on the radii of M- 
dwarfs is therefore a key component in our understanding 
of these stars. Furthermore, it is a limiting factor in char- 
acterising the planetary companions of M-dwarfs, which are 
arguably our best targets in the search for habitable worlds 
and the study of other Earth-like atmospheres. 

More measurements of the masses, radii and orbital pe- 
riods of M-dwarf eclipsing binaries, spanning both the fully 
convective regime and partially convective mass regime, for 
active and non-active stars, across a range of periods ex- 
tending beyond 5 days, are necessary to provide stringent 
observational constraints on the role of activity in the evo- 
lution of single low-mass stars. However, the influence of 
spots on the accuracy to which we can determine the radii 
from light curves will continue to impe de these efforts, even 
in the most care ful of cases (see e.g. IMorales et alj 120101 ; 
llrwin et al.ll201ll ). 

This work has studied only one third of the M-dwarfs 
in the WFCAM Transit Survey. Observations are on-going 
and we expect our catalogue of M-dwarf eclipsing binaries 
to increase. This forms part of the legacy of the WTS and 
will provide the low-mass star community with high-quality 
MEB light curves. Furthermore, the longer the WTS runs, 
the more sensitive we become to valuable long-period, low- 
mass eclipsing binaries. These contributions plus other M- 
dwarf surveys, such as MEarth and PTF/M-dwarfs, will ul- 
timately provide the observational calibration needed to an- 
chor the theory of low-mass stellar evolution. 
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Name 



HJD 



JWTS 
(mag) 



(mag) 



19a-l-02980 2454317.82863842 14.846761 0.004826 
19a-l-02980 2454317.84010834 14.844511 0.004894 



Table A2. The WTS J-band light curves for the remainder of the 
WTS MEB catalogue given in Table [All Magnitudes are given in 



for other systems. The errors, crj, are estimated using a stan- 
dard noise model, including contributions from Poisson noise in 
the stellar counts, sky noise, readout noise and errors in the sky 
background estimation. (This table is published in full in the on- 
line journal and is shown partially here for guidance regarding its 
form and content.) 

APPENDIX A: THE FULL WTS 19HR FIELD 
M-DWARF ECLIPSING BINARY SAMPLE 

In Table |A"T1 we present the periods, epochs, effective tem- 
peratures, J-band and i-band magnitudes of the 13 remain- 
ing 19hr detached, well-sampled M-dwarf eclipsing binaries 
found with this study (J ^ 16 mag). The temperatures are 
based on the SED fitting described in Section l4~T1 and may be 
under-estimated. The periods and epochs are based only on 
least-square fitting which under-estimates the errors. These 
results are accurate to ~ 30 minutes and we recommend 
to anyone planning to observe these objects in a time crit- 
ical manner that they check these values themselves with 
the light curve data provided with this paper. Note that 
19g-4-02069 is the subject of a near future publication (Nefs 
et al. in prep.) using RVs follow-up already obtained with 
GNIRS / GEMINI. The phase-folded light curves are shown 
in Figures [All and lA2l and the light curve data are provide 
in Table El 

Table IA3I contains the literature data used to create 
Figures [Toll 111 and 1121 The literature data was selected with 
the following filters: mass errors < 6.4% and radius errors < 
5.5% (comparable to or better than the errors we presented 
for the three characterised MEBs in this paper), and in the 
range 0.19 < K'h 0.71 and 0.19 < R* «J 0.71. 
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Name 


RA 


Dec A^epocha RMS 


P 


To 


J (Vega) 


i (Vega) 


T e ff,SED 




(deg) 


(deg) (mag) 


(days) 


(HJD) 


(mag) 


(mag) 


(K) 



19a-l-02980 


292 


71276 


36 


312725 


893 


5 


8 


2 


103525 


2454318 


65422 


11 


861 


± 





004 


16 


166 


± 





004 


3946 


± 


100 


19c-3-08647 


294 


30659 


36 


815037 


893 


15 








867466 


2454318 


50614 


14 


812 


± 





004 


16 


171 


± 





004 


3883 


± 


100 


19c-4-11480 


293 


81149 


36 


902880 


893 


20 


4 





681810 


2454317 


89071 


15 


850 


± 





006 


17 


208 


± 





007 


3946 


± 


100 


19d-2-07671 


294 


58622 


36 


386467 


891 


48 


9 





614540 


2454317 


99692 


15 


971 


± 





007 


17 


101 


± 





007 


4209 


± 


100 


19d-2-09173 


294 


50246 


36 


365239 


891 


22 


1 


3 


345469 


2454320 


15668 


15 


185 


± 





005 


16 


343 


± 





005 


4209 


± 


100 


19e-2-02883 


293 


32813 


36 


241312 


898 


10 


6 





810219 


2454317 


90290 


15 


976 


± 





007 


17 


272 


± 





007 


3946 


± 


100 


19f-l-07389 


292 


89403 


36 


143865 


904 


18 


3 





269868 


2454317 


97411 


15 


504 


± 





005 


16 


575 


± 





005 


4209 


± 


100 


19f-4-05194 


292 


81253 


36 


590539 


904 


35 








589530 


2454318 


10730 


16 


013 


± 





007 


17 


070 


± 





006 


4209 


± 


100 


19g-l-13215 


293 


63655 


36 


249009 


898 


10 


2 


2 


843515 


2454318 


34495 


15 


985 


± 





007 


17 


589 


± 





008 


3374 


± 


100 


19g-2-08064 


294 


16931 


36 


162723 


898 


14 


8 


1 


720410 


2454317 


94781 


11 


466 


± 





003 


15 


596 


± 





004 


4209 


± 


100 


19g-4-02069 


293 


76480 


36 


521247 


898 


11 


2 


2 


441759 


2454321 


78532 


14 


843 


± 





004 


16 


911 


± 





006 


3054 


± 


100 


19h-2-00357 


294 


66466 


36 


272874 


885 


8 


3 


7 


004082 


2454320 


79766 


15 


531 


± 





005 


16 


808 


± 





006 


3946 


± 


100 


19h-2-01090 


294 


62103 


36 


262345 


886 


11 


5 


5 


285051 


2454322 


78131 


15 


681 


± 





006 


16 


843 


± 





006 


4209 


± 


100 



Table Al. The first release of the WTS M-dwarf Eclipsing Binary Catalogue detailing the remaining MEBs in the WTS 19hr field with 
J ^ 16 mag that are not characterised in this paper. Note that 19g-4-02069 is the subject for a near future publication by Nefs et al. (in 
prep.) using RV follow-up from GNIRS/GEMINI. Please see appendix text for caveats on the quoted ephemerides. 
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Figure Al. Phase- folded light curves of the MEBs discovered in the WTS 19hr field with J ^ 16 mag... 
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Figure A2. cont... Phase-folded light curves of the MEBs discovered in the WTS 19hr field with J ^ 16 mag. 
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Name 


Period 


Mass 




Radius 


(7 r> 


-*- (..'11 


CTT cff 


Ref 




f'da.vs'l 


(Mn) 


(Mr.) 


(Rn) 


(Rn) 


(K) 


(K) 




MEBs 


















NSVS01031772A 


0.368 


0.5428 


0.0028 


0.5260 


0.0028 


3614.1 


67.2 


(1) 


NSVS01031772B 


0.368 


0.4982 


0.0025 


0.5087 


0.0031 


3515.6 


32.5 


(1) 


GUBooA 


0.489 


0.6100 


0.0071 


0.6230 


0.0163 


3917.4 


128.3 


(1) 


GUBooB 


0.489 


0.5990 


0.0061 


0.6200 


0.0203 


3810.7 


133.9 


(1) 


MG1-1819499A 


0.6303135 


0.557 


0.001 


0.569 


0.002 


3690.0 


100.0 


(2) 


MG1-1819499B 


0.6303135 


0.535 


0.001 


0.500 


0.003 


3610.0 


100.0 


(2) 


GJ3236A 


0.77126 


0.376 


0.016 


0.3795 


0.0084 


3312.0 


110.0 


(3) 


GJ3236B 


0.77126 


0.281 


0.015 


0.300 


0.015 


3242.0 


108.0 


(3) 


YYGemA 


0.814 


0.5974 


0.0047 


0.6196 


0.0057 


3819.4 


98.0 


(1) 


YYGemB 


0.814 


0.6009 


0.0047 


0.6035 


0.0057 


3819.4 


98.0 


(1) 


MG1-116309A 


0.8271425 


0.567 


0.002 


0.552 


0.004 


3917.4 


100.5 


(2) 


MG1-116309B 


0.8271425 


0.532 


0.002 


0.532 


0.004 


3810.7 


97.8 


(2) 


CMDraA 


1.268 


0.2310 


0.0009 


0.2534 


0.0019 


3133.3 


73.0 


(1) 


CMDraB 


1.268 


0.2141 


0.0009 


0.2396 


0.0015 


3118.9 


102.2 


(1) 


MG1-506664A 


1.5484492 


0.584 


0.002 


0.560 


0.001 


3732.5 


104.6 


(2) 


MG1-506664B 


1.5484492 


0.544 


0.002 


0.513 


0.001 


3614.1 


101.3 


(2) 


MG1-78457A 


1.5862046 


0.5270 


0.0019 


0.505 


0.008 


3326.6 


101.1 


(2) 


MG1-78457B 


1.5862046 


0.491 


0.002 


0.471 


0.009 


3273.4 


99.5 


(2) 


LP133-373A 


1.6279866 


0.34 


0.02 


0.330 


0.014 


3144.0 


206.0 


(1) 


LP133-373B 


1.6279866 


0.34 


0.02 


0.330 


0.014 


3058.0 


195.0 


(1) 


MG1-646680A 


1.6375302 


0.499 


0.002 


0.457 


0.006 


3732.5 


51.9 


(2) 


MG1-646680B 


1.6375302 


0.443 


0.002 


0.427 


0.006 


3630.8 


50.5 


(2) 


MG1-2056316A 


1.7228208 


0.4690 


0.0021 


0.441 


0.002 


3459.4 


179.8 


(2) 


MG1-2056316B 


1.7228208 


0.382 


0.002 


0.374 


0.002 


3318.9 


172.5 


(2) 


KOI126B 


1.76713 


0.2413 


0.0030 


0.2543 


0.0014 


- 


- 


(5) 


KOI126C 


1.76713 


0.2127 


0.0026 


0.2318 


0.0013 


- 


- 


(5) 


HIP96515Aa 


2.3456 


0.59 


0.03 


0.64 


0.01 


3724.0 


154.0 


(0) 


HlP9o515Ab 


2.3456 


0.54 


0.03 


0.55 


0.03 


3589.0 


157.0 


(6) 


CUCncA 


2.771 


0.4333 


0.0017 


0.4317 


0.0052 


3162.3 


156.7 


(1) 


CUCncB 


2.771 


0.3980 


0.0014 


0.3908 


0.0095 


3126.1 


154.9 


(1) 


1RXSJ154727A 


3.5500184 


0.2576 


0.0085 


0.2895 


0.0068 






(7) 


1RXSJ154727B 


3.5500184 


0.2585 


0.0080 


0.2895 


0.0068 


- 


- 


(7) 


LSPMJ1112A 


41.03236 


0.3946 


0.0023 


0.3860 


0.005 


3061.0 


162.0 


(8) 


LSPMJ1112B 


41.03236 


0.2745 


0.0012 


0.2978 


0.005 


2952.0 


163.0 


(8) 


Kepler 16 A 


41.079220 


0.6897 


0.0035 


0.6489 


0.0013 


4450 


150 


(9) 


Keplerl6B 


41.079220 


0.20255 


0.00066 


0.22623 


0.00059 






(9) 


Non-M-dwarf primary EBs 


















NGC-2204-S892B 


0.4520000 


0.6621 


0.0050 


0.6800 


0.0203 


3944.6 


110.5 


(1) 


IM-VirB 


1.3090000 


0.6644 


0.0048 


0.6809 


0.0131 


4246.2 


129.0 


(1) 


RXJ0239B 


2.0720160 


0.693 


0.006 


0.703 


0.002 


4275.0 


109.0 


(10) 


MD-WD EBs 


















RXJ2130 


0.5210356 


0.555 


0.023 


0.553 


0.017 


3200.0 


100.0 


(10) 


Interferometry 


















GJ411 




0.403 


0.020 


0.393 


0.008 


3570.0 


42.0 


(11) 


GJ380 




0.670 


0.033 


0.605 


0.020 






(11) 


GJ887 




0.503 


0.025 


0.459 


0.011 


3797.0 


45.0 


(12) 



Table A3. Literature values for systems used in Figures 1101 1111 and 1121 with mass errors < 6.4% and radius errors < 5.5%, in the 
range 0.19 ^ ^ 0.71 and 0.19 ^ -R* ^ 0.71. Temperatures are given when available in the literature, but those without are not 
included in Figure 1101 There are no rotation periods given for the interferometric measurements therefore these ar e ex cluded from 
Figure IT2l References: (1) DEBC at a nd references therei n (w ww .astro . keele . a c . u k/.jkt/de bcat/l, (2llKraus et al] <201lD. (31 1 Irwin et aU 
[|2009ll, (41 IVaccaro et ail ll2007f). (51 ICarter et al] feOll]) , (6) iHuelamo etail d2009ll, (71 lHartman et alj <201lll. (81 llrwin et al.l <201lll . 
Ol lDovle et al] j201 llT TlOl iKnigge et al] l|201ll l and references therein. (Ill ISegransan et al] l|2003r i. (12l lDemorv et al] l|2009l l. 
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